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CALCULATED SPANWISE LIFT DISTRIBUTIONS AND AERODYNAMIC 
INFLUENCE COEFFICIENTS FOR UNSWEPT WINGS 
IN SUBSONIC FLOW 


By Franklin W. Diederich and Martin Zlotnick 
SUMMARY 


Spanwise Lift distributions have been calculated for nineteen 
unswept Wings with various aspect ratios and taper ratios and with a 
variety of angle-of-attack or twist distributions, including flap and 
aileron deflections, by means of the Weissinger method with eight control 
points on the semispan. Also calculated were aerodynamic influence coef- 
ficients which pertain to a certain definite set of stations along the 
span, and several methods are presented for calculating aerodynamic 
influence coefficients for stations other than those stipulated. 


The information presented herein can be used in the analysis of 
untwisted wings or wings with known twist distributions, as well as in 
aeroelastic calculations involving initially unknown twist distributions. 


INTRODUCTION 


In the design and development of an airplane, a knowledge of the 

spanwise lift distribution on the wing is important in predicting the 
structural loads and the stability characteristics. For high-speed air- 
planes having flexible wings, the calculation of the spanwise lift dis- 
tribution is an aeroelastic rather than a purely aerodynamic problem. 
In aeroelastic calculations means are required for calculating the span- 
Wise lift distribution for angle-of-attack (or twist) distributions which 
are initially unknown. Aerodynamic influence coefficients constitute the 
most convenient of these means. 


One of the most satisfactory techniques developed in recent years 
for calculating the spanwise lift distribution on a wing in subsonic flow 
has been the Weissinger L-method (ref. 1), which can be applied to a 
large variety of plan forms and yields solutions of sufficient accuracy 
for all practical purposes without requiring an unduly long time for the 
calculations. This method may be considered as a simplified lifting- 
surface theory because the calculation of the lift on the wing is treated 
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as a boundary-value problem, the boundary condition being that the down- 
wash angle induced by the bound and trailing vortices is equal to the 
geometric angle of attack at the three-quarter-chord line. 


In the present paper, symmetrical and antisymmetrical lift distri- 
butions and some associated aerodynamic parameters have been calculated 
by means of the Weissinger method with eight control points on the semi- 
span for several continuous and discontinuous angle-of-attack conditions 
on nineteen unswept wings having various aspect-ratios and taper ratios. 
A convenient matrix formulation of the Weissinger method was used in con- 
junction with the Bell Telephone Laboratories X-66744 relay computer at 
the Langley Laboratory to make the calculations. This formulation is 
described in appendix A. 


Aerodynamic influence coefficients have been calculated for these 
nineteen wings for a certain prescribed set of stations along the span 
and are presented herein, and several methods for calculating aero- 
dynamic influence coefficients for any arbitrary set of stations from 
the numerical results of this paper are also presented. The influence 
coefficients calculated in this manner can be used in aeroelastic anal- 
yses similar to that of reference 2. 


SYMBOLS 
A aspect ratio 
b wing span ` 
bail aileron span 
be flap span 
Com root bending-moment coefficient for unit angle of attack, 
2 x Bending moment 

qSb 
Cp, induced-drag coefficient at a unit angle of attack 
Cr, Lift coefficient at a unit angle of attack 
στο lift-curve slope per radian for additional-type loading 
τη ур lift coefficient for half of antisymmetrically loaded wing at 


1/9 

a unit tip angle of attack, =g 
q = 
2 


L. 
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Rolling moment 


C1 rolling-moment coefficient, 
qSb 
Clg = E 
“Ly coefficient of damping in roll 
Ств/ 05 rolling-moment coefficient for unit aileron deflection 
с wing chord 
c average chord, S/b 
е; section lift coefficient 


integrating matrix for Сам (see appendix A) 
integrating matrix for Сү (see appendix А) 
integrating matrix for бтлу (see appendix А) 


integrating matrix for С) (see appendix A) 


lift on one semispan 


free-stream Mach number 
aerodynaemic-influence-coefficient matrix 
dynamic pressure 

wing area 

free-stream velocity 

lateral ordinates 


spanwise location of discontinuity 


spanwise center-of-pressure location 


angle of attack, radians 
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de, /45 
ας effective angle of attack for unit flap deflection, ------ 

de, /до 
r vortex strength . | τ | і жш 
p dimensionless vortex strength, z NS сХс; 

bV 
Š flap or aileron deflection angle, radians 
Ө = cost y* 
ER -1 x 

8, = cos Yo 


A sweepback angle, deg 
À taper ratio 
Subscripts: 

a antisymmetrical 

ail aileron 

C continuous 

D discontinuous 

f flap 

L left | | | τος 
Β right 

S symmetrical 

t tip (7 
Superscript: 

Е 


dimensionless with respect to semispan »b/2 
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Matrix notation: 


| | row matrix 
1 | column matrix 
| | general matrix (not a row or a column matrix, but need not 


be square ) 
| | diagonal matrix 


BP unit (identity) matrix 
In matrix notation, a prime indicates the transpose of the matrix. 


PRESENTATION OF CALCULATED RESULTS 
Spanwise Lift Distributions 
Geometric characteristics. ОҒ the nineteen plan forms treated in 
this paper are indicated in table I. Lift distributions due to the 


following continuous symmetric and antisymmetric angle-of-attack dis- 
tributions have been calculated for each of these plan forms: 


Symmetric angle-of-attack distributions: 
Constant (o = 1) 
Linear (0 = |у%|) 
Quadratic (о = y*?) 
Cubic С = | #3 |) 
Straight-line (а - уж) 
Antisymmetric angle-of-attack distributions: 
Linear (о = y*) 
Quadratic (a = үке for y* 20; а = -уж2 for ух% °) 
Cubic (o = y*? ) 
Quartic (a = уж for y* 20; a= yx} for y* < o) 
Quintie (а = уж?) 
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The straight-line angle-of=attack condition was included to repre- 
sent actual structural twists where the surface of the wing is generated 
by straight lines so that the product c*a, the deflection of the leading 
edge, varies linearly with y*; that is, 


c*a = ety“ 


or, for unit twist at the tip, 


For untapered wings, the straight-line lift distribution is the same as 
the linear lift distribution, and for wings of zero taper ratio, it is 
undefined. | 


Lift distributions for flap-type and aileron-type angle-of-attack 
distributions are also presented. A correction which has been made for 
the spanwise discontinuity in the angle of attack is derived in 
appendix B. Тһе values of— bp/b and b,,;/b (ratios of the flap span 


to the total span and the aileron span to the total span, respectively ) 
for which the lift distributions have been calculated are 0.1, 0.2, 0.3, 
0.1, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.0. As is usual, the flaps have been 
taken to be inboard and the ailerons outboard. The lift distribution 
for any flap or aileron configuration may be obtained, however, by 
linear superposition; thus, the lift distribution for an outboard flap 
extending, for example, from y* = 0.5 to y* = 1.0 can be obtained 


b 
by subtracting the lift distribution for the inboard flap E. o.5) 


b 
from the additional lift distribution (= 5 1.0). A similar procedure 


can be used for inboard ailerons. 


The lift distributions pertaining to each plan form are given in 
one figure, parts (&) and (b) showing the lift distribution due to sym- 
metric and antisymmetric continuous angle-of-attack distributions, - 
respectively, and parts (c) and (d) showing the lift distribution due 
to flaps and ailerons, respectively. Lift distributions for plan forms 
with am aspect ratio approaching zero have been taken from reference 3 
and are included herein in figure-1 for the sake of completeness. (As 
indicated in ref. 5, the lift distribution for a wing of very low aspect 
r&tio is independent of the plan form, provided the trailing edge is not 
reentrant.) The lift distributions on the nineteen wings considered in 
the present p&per &re presented in figures 2 to 20. Table I serves &s 
а table of contents for this group of figures. 
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Aerodynamic Parameters 


The aerodynamic parameters ы! Сам; у; бр. C14: and. "I4 уо 


for the nineteen plan forms considered are compiled in table ІІ. Тһе 
values of Ст and Срм for a unit effective flap deflection are pre- 


sented in table III, and the values of στη уа and C15 /% for a unit 


effective aileron deflection are presented in table IV. 


Aerodynamic Influence Coefficients 
for Stipulated Stations 


Aerodynamic influence coefficients for symmetric and antisymmetric 
lift distributions were obtained as shown in appendix A and are pre- 


sented as the matrices [Še | and [д | in tables V(&) and V(b), respec- 


tively. Each influence-coefficient matrix in the table applies to a 
given plan form. These influence-coefficient matrices are used to cal- 
culate the spanwise lift distribution for any continuous angle-of-attack 
condition from the following matrix expressions: 


{Pep = Cry [8] {ов} 
{Pay = Cra [йы] {oa} 


for the symmetrical and antisymmetrical distributions, respectively, 
where α is the angle of attack at stations y* = 0.9808, 0.9250, 
0.8515, 0.7071, 0.5556, 0.3827, 0.1951, and О and ΓΈ is the desired 
lift at these stations. Іп this paper the convention is that the angle 
of attack for the station nearest the wing tip (y* = 0.9808) is the 


first element of the angle-of-attack matrix {a} and the lift at the 
same station is the first element of the lift-distribution matrix {γ᾽}, 
The matrices [as] and [aa] of table V are arranged accordingly. 


DEVELOPMENT OF AERODYNAMIC INFLUENCE COEFFICIENTS 


FOR ARBITRARY STATIONS 


Ihe influence coefficients described in the preceding section are 
satisfactory for many purposes; for instance, the stipulated stations at 


8 NACA TN 3014 


which the lift is given are convenient for plotting spanwise lift distri- 
butions because the points are concentrated near the wing tip where the 
curvature of the lift distributions is greatest. In some cases, however, 
other considerations may determine the points on the span at which the 
lift is to be calculated. For instance, when the influence coefficients 
are to be used in an aeroelastic analysis, the location of the stations 
may be dictated by the structural characteristics of the wing; also, if 
lift distributions are to be calculated for the sake of comparison with 
experimental results, this comparison can be facilitated by c&lculating 
the lift at the same stations at which it is measured and thus avoiding 
the necessity of graphical or numerical interpolation. 


The following sections describe several methods for developing 
aerodynamice influence coefficients for arbitrary stations. 


Method Using Interpolating Matrices 


One way of constructing an influence-coefficient matrix for any 
stations from the matrices presented herein is to calculate interpo- 
lating matrices which give the angles of attack at y* = 0.9808, 0.9239, 
0.8515, 0.7071, 0.5556, 0.3827, 0.1951, and О in terms of the angles of 
attack at the given stations and the values of the lift at the given 
stations in terms of those &t the stations y* = 0.9808, 0.9259, 0.8315, 
0.7071, 0.5556, 0.5827, 0.1951, and O. Тһе desired influence-coefficient 
matrix would then be obtained by postmultiplying the one given herein by 
the angle-of-attack interpolating matrix and premultiplying it by the 
Lift interpolating matrix. 


In order to illustrate the nature of these calculations, a pair of 
matrices are calculated for stations at every tenth of the semispan by 
linear interpolation. (The stations at-which the angle of attack is 
given and those at which the lift is to be found need not be the same, 
but here, as in most cases, they are chosen to be the same as a matter 
of convenience. ) 


With linear interpolation, 
20.9808 = 0.9080 0% 0.19200, ο 


00.8515 = 0.31509, ο + 0.68500. в 


Up = Sp d 
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which can be written in matrix form as 


a αρ. 9808 0.808 0.192 ο 0 0 0 
αρ, 9239 0.259 0.761 ο 0 0 ο 
αρ 8515 O 0.315 0.68 ο ο ο 
αρ, 7071 ο O 0.071 0.929 ο ο 
09.5556 ο o o о 0.556 о. 
00.5827 0 0 0 0 0 о 
αρ 1951 0 0 0 0 0 0 
оо 0 0 0 0 0 0 


0 


ο 


0.827 0.175 


0 


0 


0 


ο 


ο 


0 


о 


0.951 0.09 


0 


о 


0 


where the rectangular matrix on the right side of the equation is the 


desired angle-of-attack interpolating matrix. 


Similarly, for the lift distribution, 


πιο” 0 
" (This information is known from physical considerations; to calculate 
πο by extrapolation from Г^о. 9808 and T* 0.9239 would give a 
қ spurious үзіле.) 
x B ο ) 
l'o.o = P"o.8515 + 0 бек 0.9259 - Г 0.8515 


0.0929 x 
Го.8 = Го. тот + τοί ο. 8515 zb 0.7071) 


a = 0.1770. 9315 - 0.2550. 7071 


% ә 8 a 
X 
го 


ні 

ж 
© 
И 


O. Tl o, 9239 + 0.259Г70. 8515 
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or, in matrix form, 


Гб 0 ο ο ο ο 0 ο Г. 9808 
T^0.9 0.781 0.259 0 ο ο ο ο ro, 9259 
г®* 8 о 0.717 0.255 0 ο ο ο το 8515 
Το, ή ο о 0.953 0.047 ο 0 0 Гб, 7071 
r*0.6 0 O 0.293 0.707 0 0 ο г*0.5556 
Do. ο ο о 0.679 0.321 0 ο Г*0.5827 
pea 0 0 о 0.100 0.900 0 ο "9.1951 
Го. ο ο 0 о 0.559 ОЛА ο ο 
r"0.2 O 6) 0 6 0.026 0.974- O 

roi ο 0 ο ο о 0.515 0.487 

ος ο ο ο ο ο о 1.000 


where the rectangular matrix on the right side of the equation is the 
desired lift interpolating matrix. 


Although linear interpolation is by far the simplest type, the 
results obtained with it are not so accurate as a higher-order interpo- 
lation procedure. Parabolic interpolation should be satisfactory for 
the angle-of-attack matrix and for most of the lift-distribution matrix, 
except near the wing tip. The numerical factors required for parabolic 
interpolation can be calculated by means of Lagrange's interpolation 
formula for an mth degree polynomial | 


m (^ - уч) 
г" (у - у 
mi | z h i 
“у%)- N (ушы 


м 
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where Ғ(у%) represents either the lift or the angle-of-attack distri- 
bution, у is the station at which f(y*) is to be determined by 


interpolation, and у, у%, + . + аге the stations at which £(y*) 


is presumed to be known. The prime mark on the product signs is the 
conventional designation of the fact that the term for 1 =k is to be 
omitted. For parabolic interpolation (m = 2) this formula reduces to 


“νὴ PET σὰ 
: (у% - Yt) - Υπ), 

MOS - Y*J(y*s - M 

(h р у% )(у%, - у%2 

MES -Yy*Sy ys, - Е 


Near the wing tip the lift distribution cannot be approximated 
accurately by an ordinary parabola but can be represented instead by & 


linear superposition of the two functions (i - yx) 2 and (1 - уж)2/2 
as suggested by V. М. Falkner. With this approximation the desired inter- 
polating factors for y* between 0.9259 and 1 are the two elements of the 
row matrix obtained by postmultiplying the row matrix by the square matrix 
in the equation 


9.652 -1.223 I". 9808 
* Z " 1/2 u 3/2 
r (у%,) (6 у%) (1 ух) να кені yan 


For y% = 0.96, for instance, the f&ctors obtained in this manner are 
0.9157 and 0.2651, so that 


г*0.96 = 0.915715. 9808 + 0-26510"0, 9239 


lagrange's general interpolation formula can be used for higher- 
order interpolation (m > 2), but the effort entailed in calculating the 
interpolating factors is not generally justified by the increase in 
accuracy obtainable compared to parabolic interpolation. 
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Interpolating factors for Oc and T* сап also be obtained by 
representing these functions by Fourier series in 8. For the lift dis- 
tribution the trigonometric interpolation formula 


г” (η) = 


ft віп jð, sin jð 
n + 1 i=l ( 1) 4-1 i h 


can be used, where n = 15 for the calculations in this paper, and 
where, аз in this paper, the values of 9 are at equal increments 


94 = 2 . For the angle-of-attack distribution this formula is not 
1 n+ 1 


applicable, because the angle-of-attack distribution cannot be repre- 
sented accurately by a finite sine series and because the angle-of- 
attack values are presumed to be given at nonequal increments іп 981. 
(If they were given at eight equal increments no interpolation would be 
required.) Although in principle a matrix-inversion method could be | 
based on an expansion of the angle-of-attack distribution in a cosine 
series, the matrices to be inverted are generally ill-behaved so that 
the results are likely to be of doubtful accuracy. 


Methods Using Green's Function 


Basis of methods.- Although the aerodynamic-influence-coeffictent 
matrices discussed in the preceding section have the property that when 


postmultiplied by the angle-of-attack matrix (a and multiplied by the 


lift-curve slope they yield the lift-distribution matrix {τ}, their 


individual elements have no direct physical significance. On the other 
hand, & structural influence coefficient has the significance that it 
represents the deformation at one point caused by a unit concentrated 
force at another point. A corresponding type of aerodynamic influence 
coefficient would represent the lift at one point y* due to a "unit 
concentrated angle of attack" at another point yg*. This angle of 


attack is actually an angle-of-attack distribution represented by a 
Dirac delta (impulse) function of the distance along the span, that is, 
a function which, in the limiting case as Ay,* approaches zero, is 


zero everywhere except in the interval y,* S ух £ у ® + Ayg* where 


the ordinate is equal to ше the area under this function would 
ο 

always be 1, which justifies the use of the term "unit" in connection 

with this distribution. The desired influence coefficients would then 

be the values (at various values of у%) of the lift-distributions due 

to these angle-of-attack distributions for various values of- yg*. 
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The desired Lift distributions, which constitute a type of Green's 
function for Weissinger's integral equation, can be obtained as follows. 
Let Г”(у%,у;%) be the lift distribution for a unit effective deflec- 


tion of a flap which is located between y* = yg* and y* = і. Then 


the lift distribution for an angle-of-attack distribution for which 
the angle of &ttack is zero everywhere except in the interval 


yg* 5 y* S yg* + Ayg*, where it is L, is given by 


Г”(ух,у%) ш D (yr ,у„®*+Лу*) 
Дус“ 


The lift distribution corresponding to the unit concentrated angle-of- 
attack distribution therefore is the limit of this expression ав Ay,* 
x 
Sr (ys yos) 
approaches zero, which is UAE о by definition. For any given 
о 
angle-of-attack distribution the lift distribution can be determined by 
linear superposition of lift distributions of the Green's function type 
ag follows: 


NM 1 Әг* * 
r*(y*) = -/ α(γοχ)------------ ay (2) 
En Sy 
The desired Green's function can thus be obtained by απο πόας. 
р *(yz,yo*) and taking its partial derivative with respect to Vo . А 


more convenient approach, however, is to consider symmetric and anti- 
symmetric loadings separately. By a repetition of the preceding argu- 
ment the following results are then obtained: 


J: a* (y* *) 
* βίου Уо 
pw (уж) = [^ (у аур, (5) 
ο ду 
1 SP” ari yo*) 
x aili Jo 
PS = - [| Aa (y 9X) —————— dy,,* (+) 
where, as elsewhere in this paper, I“, is the lift distribution (as а 
function of y.) for an inboard flap extending over the interval 
-у Sy" S E % and I". у is the lift distribution for outboard 
ailerons extending over the intervals yo“ < |у*| < 1. The argument 


ус in a(yg*) can be regarded simply as a variable of integration 
corresponding to у^, 
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The desired Green's function can thus be obtained from the flap and 
aileron distributions given herein by differentiation with respect io 
Уо*. The results presented in parts (c) and (d) of figures 1 to 20 can, 
for instance, be cross-plotted as functions of γον (νο» ==: for gym- 
Pail 

b 
given values of y* and the differentiation then performed graphically. 
The graphical procedure, however, is tedious and relatively inaccurate. 
Similarly, numerical differentiation of these lift distributions would 
tend to be inaccurate. Two numerical procedures which avoid differenti- 
ation of the lift distributions are therefore presented in the following 
sections; one consists in calculating the desired Green's functions 


directly, and the other consists in using derivatives of the angle-of- 
attack distribution. 


metrical loadings, and y. = 1 - 


for antisymmetrical loadings) for 


Direct calculation of Green's functions.- Inasmuch as the desired 
Green's functions are lift distributions corresponding to angle-of- 
attack distributions defined by delta functions, they can be calculated 
directly provided the singularities in the engle-of-&ttack and lift dis- 
tributions are taken into &ccount. Appendix B of this p&per describes 
the method by which the singularities in the flap- and aileron-type 
angle-of-attack and lift distributions were taken into account. This 
method is extended to the case of impulse-type angle-of-attack distri- 
butions in appendix C. 


The resulting lift distributions I*,' and Γρ are identical 


or Yt y у%,у 
with the Green's functions ЙІ; ы/а and бауы) They 


Зуо* дуо* 
contain logarithmic singularities at у.* = y* which must be split off 


before the integrals in equations (5) and (4) can be evaluated numeri- 
cally. For a symmetric distribution the procedure is as follows. For 
a given value of у*, 


2. 
r* (y*) = J. [as(vo*) - eg Gr) |r*s (уу ау қ + 


1 1 
σρίν») | rhe (tore) 6) 


The first ая does not contain any singularity; the integrand is 
zero at y. = y*. Тһе second integral canm be evaluated explicitly and 


is 
ας OT. (у%) 
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where Icy) is the lift distribution for а = 1 over the entire 


span, so that the right side of equation (5) cam be evaluated numerically 
without difficulty by using any set of integrating factors appropriate to 
the stations of interest, such as those of Simpson's rule if the points 
are equally spaced and the number of intervals between y* = O and 

ух = 1 is even. 


If these integrating factors are written in the form of a diagonal 
matrix and designated by |1| , equation (5) can be written аз 


ге G) = [Te (res) |l {οα(νον) - о.б) + αρίνκ)γ᾽ε, ὧν’) (6) 
Now, if the row matrix | Lyx | is defined as 


|з|-|0 0000100... οἱ 


with the element 1 at the position corresponding to y* and zeros 
elsewhere, if the matrix E is defined as & square matrix all the 


rows of which are equal to | ty* |; and if (11 represents the unit 


matrix, then 


Es Gen] un er - Be] | ne 
r*, қ (у%) Ем {ας (ус®*) 


СЕТИ [rs (уус) LE 5 | Ей! {бе} 


г“ (уж) 


or 


г (y*) = ств] , 19875) (7) 


where the matrix ЕЛЕР defined by 


[| x _ тыге ТЕ г” gen Б И - EMI (8) 
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is the row matrix corresponding to y* of the desired influence- 
x 
C 


directly from the curves labeled "Constant" of parts (a) or the curves 


coefficient matrix (941. The v&lues of . сап be obtained 


b 
for Te 1.0 of parts (c) of figures 1 to 20, and the calculation of 
ж 
lg (yt y e) 
“Ley 
repeated for all other values of уў to obtain ба, so that, finally 


uan (y*) r*a (у%у%) 
[as] -|-x— |+ EC wm LN (9) 


where P | is а diagonal matrix in which the elements аге the sums of 
8 


is described in appendix C. This calculation must be 


the elements in the rows of the matrix 
X 1 
r*a (y*,yo*) үг] 
"Ig 


Similarly, for an antisymmetric distribution the singularity can be 
split off in several ways, one of them being the following. For & given 
value of .y* equation (4) can be written as 


(сы 1 : M d qd "S | * (уж) 
Pam) =f” [ρο") - аг, Grot νο” МАСА 


where a (у*) is the lift distribution for & unit effective displace- 


ment ОҒ a full-span aileron. Іп matrix notation this relation may be 
written as 


га) = Са, [84]. , Jos (уо*) (10) 


3A 
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where 
Е тен мб тым) ) [τ] НЕ (11) 
[Sal 1-6 сы ЕИ 
The values of о can be obtained from the curves for —— = 1.0 
la 


of parts (d) of figures 1 to 20. This procedure must again be repeated 
for each value of y* to obtain all the rows of [Q,|, with the result 


that 
RJ- - г ue : борий | 5 а 


where the diagonal matrix р | consists of the sums of the elements 


of the rows of the matrix 
Li 
Сі 


ΙΤ] 


d 


In order to indicate the extent of the effort involved in calcu- 
lating these influence-coefficient matrices, a step-by-step summary of 


these calculations is given for | 9% ; with the obvious modifications 
this procedure also applies to | 3α|- In the first six steps Green's 


functions are calculated in accordance with the procedure indicated in 
appendix C; in the remaining steps | | is calculated in accordance 
with the previous discussion in this section. 

(1) For the values of 0. = сов-1 y,* of interest, the values of 
sin обо (for n=1, 3,5, + + . 15) are obtained from trigonometric 
tables and assembled in a matrix 4sin n9 o 
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-L 
(2) This matrix is premultiplied by the matrix | [eos ΠΝ : 
(Зее table VI.) 


(5) The resulting matrix is premultiplied by the matrix ЕЛ con- 


taining the elements  F(n*,y*) defined in appendix A. The evaluation 
of this [Fe | matrix is probably the most time-consuming part of the 


calculation because it does not lend itself very readily to high-speed 
automatic computation. For any one of the nineteen plan forms considered 


in this paper, the matrices ағы and 2, | are available upon request 
from the Nationsl Advisory Committee for Aeronautics. 


(4) The resulting matrix is premultiplied by the matrix [ας | given 


in the present paper for the plan form of concern, and the matrix obtained 
in this manner is multiplied by the constant 1/π. 


(5) The values of ea are calculated from equation (C3) for the 


given values of Ө. and for 9 "ES n = 1, 2, . .. 8; they are then 


divided by τ. and assembled in a matrix the columns of which 


pertain to given values of 9,5. 


(6) The matrices obtained in steps (4) and (5) are added to obtain 


ж? 
the matrix E 
Ly 
* 
r gi (ΥΣ) 
(7) Values of cw are obtained by reading the values of 
ес, 
— at y* defined by y* = cos nA n = 1, 2, . . . 8, from figures 2 
to 20 for a constant value of α and multiplying them by =, these 


2 
Us. (y*) 
Cr, | 


values аге assembled in а diagonal matrix 
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(8) The matrix obtained in step (6) is postmultiplied by the 
diagonal matrix of integrating factors B3? 


(9) The di&gonal matrix > | is assembled from elements calcu- 


4. adding the elements іп a given row of the matrix obtained in 
step (3). 


(10) The matrices obtained in steps (7) and (8) are added to each 
other and th&t obtained in step (9) is subtracted from the sum. The 


resulting matrix is the desired matrix las]; as defined in equation (9). 


The entire calculation thus involves the calculation of 8P values 


sin nO, 


of с апа of Га (P being the number of stations Уо*); 


four matrix multiplications, and three matrix additions, as well as the 
calculation of the 64 elements of the matrix [ғ], if this matrix is 


not available. 
Calculation using spanwise derivative of the le-of-attack distri- 
bution.- Equation (5) can be integrated by parts to yield 
l x 
x< y 
P'.(y*) = ag Qr" e (y, 1) “4 Balto" Jen (ys уа (15) 
O ο 


where I"e.(y*,1) is the lift distribution for a unit effective deflec- 


tion of a full-span flap. Similarly, integrating equation (4) by parts 
yields 


ra (y*) - Г Lm Sis ro) (у ayo σος (1%) 


Here again the argument va" in о(уо%) is merely а variable of inte- 


gration, corresponding to y*. In these equations the lift distri- 
butions τὸ г and ps ail serve as influence functions, so that no Green's 


functions need be calculated. Furthermore, neither these influence 
да(у0%) 
functions nor, in most cases of interest, the functions = have 
ο 
singularities in the range of integration, so that the numerical evalu- 
ation of the integrals of equations (13) and (14) can be effected very 
readily. 
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With a set of integrating factors [т] for the stations of interest 


E Sog (yo *) m 
and the identity a,(1) = a,(0) + E dy,*, equation (15) сап be 
О ο 


written as 


Us ont - ов (ο) r*e (9,1) E + [Em - 
REGEM: | ιτ] faen (15) 


where [τὴς (9,1) | is a matrix all the columns of which are equal to the 


column matrix ЫСЫ ‚ and | I | is the diagonal matrix consisting 


of the integrating factors. Equation (15) can be rewritten in terms of 
& new influence-coefficient matrix ΓΙ defined by 


lI 


Г г(у%у0”) [т] (16) 


(η. || Ге) | Pe 
pg = [Z]. m 


as 


{το} = os(o) ежа) + ANE (17) 


_ 
Similarly, with а new influence-coefficient matrix |6) a | defined 
by : vs š — ‘ = 


X 
ΕΝ - d ш 


ЕН (18) 
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equation (14) can be written as 


Ы (у% y - ΠΕ» (19) 


The matrices СЯ апа СЯ are based on the &ssumption that 
x 
0 (ус ) 
x 
ду 
lated when Оң is discontinuous. Discontinuities in aq result from 


is nonsingular and continuous. However, this assumption is vio- 


control-surface deflection or from deflection of parts of the wing rela- 
tive to the rest of the wing and can be treated in the manner indicated 
in appendix B or, more simply, by superposition of the Lift distributions 
given in figures 1 to 20. Іп the angle-of-attack distributions for which 
influence coefficients are particularly useful, namely those due to 
structural deformations, discontinuities cannot occur. 


ἀα(γοχ) 
дус% 
for wings with discontinuous stiffness distributions. Actually this 
theory is inapplicable for such wings, and the spanwise slope of the 
twist is never discontinuous. If, however, simple-beam theory is to be 
used anyway for the sake of convenience and because the errors involved 


are considered to be acceptable, then the matrices [a's] and Га". | 


Discontinuities in can &rise if simple beam theory is used 


can still be used provided one of the stations is located at tbe point 
of the discontinuity in the stiffness distribution and provided suitable 
integrating factors are used. 


Ihe objection may be raised against the influence-coefficient mat- 
rices [a's] and СЯ that they do not actually express the lift 


distribution in terms of the angle-of-attack distribution but rather 
require its derivative.  Inasmuch as the angle-of-attack distribution 
can always be reduced to a continuous one by splitting off the discon- 
tinuous part and treating it as described elsewhere in this report, the 
derivative of the angle-of-attack distribution can be obtained numeri- 
cally by using numerical differentiating factors obtained from any text 
on numerical analysis. These differentiating factors can be assembled 
into a differentiating matrix which when postmultiplied by the matrix 
of the angle-of-attack values yields a matrix of values of the spanwise 


derivative of the angle-of-attack distribution. The matrices ЕРІ 


and [a'a] can then be postmultiplied by this differentiating matrix 
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in order to obtain new influence-coefficient matrices which express the 
lift distribution directly in terms of the angle-of-attack distribution. 
However, the main advantage of the method outlined in this section is 
that in aeroelastic calculations, for which aerodynamic influence coef- 
ficients are primarily intended, the angle-of-attack distribution usually 
is obtained by integrating its derivative; the use of the derivative then 
actually saves a calculation. 


In such aeroelastic calculations the lift distributions can be con- 
sidered to consist of a known "rigid wing" part (due to airplane attitude 
or motion, built-in twist, or control deflection), which can be calculated 
initially with due regard to all discontinuities, and an initially unknown 


part due to structural deformation; the matrices СА and [a'a] can 


be used. to advantage in calculating the latter part. The calculation of 
a(y*) can then be obviated altogether, because if the structural defor- 
mations are referred to the plane of symmetry the structural part of 
a(y*) is zero for y* = O, so that the first term on the right sides of 
equations (15) and (17) represents a known rigid-wing lift distribution 
and can be included with the others. Thus, in general, 


τα} = {τ᾽} ici wing + cy, [8'] А а, | 


The separate treatment of these two parts in an aeroelastic analysis 
presents no difficulties and can be effected in a manner similar to that 
employed in reference 4 for the lift distribution due to aileron deflec- 


tion; the use of [@'| rather than [Q] requires only the omission of 
one of the integrating matrices in the methods of references 2 and і. 


DISCUSSION 


General Limitations of the Results of This Paper 


The Weissinger L-method, its range, and validity have been dis- 
cussed in several previous papers (for example, refs. ^ and 5) so that 
in this section only & few comments sre made about special applications. 


Number of control points.- In references l and 5, four control 
points on the semispan were found to give satisfactory accuracy for the 
lift distributions due to constant angle of attack. in the present 
paper, interest is centered primarily on the lift distributions due to 


МАСА TN 5015 25 


twist &nd control deflection, and for these cases the additional effort 
entailed in using eight rather than four control points was believed to 
be warranted by the resulting increase in accuracy. 


Fuselage, nacelle, and tip-tank interference.- The Weissinger 
method and all results presented apply only to wings without fuselages, 


nacelles, or tip tanks. At low angles of attack the lift distribution 
on the wing is not affected to a large extent by the presence of the 
fuselage except when it covers a large part of the wing; the effect is 
largely localized near the wing root and is most pronounced for the 
constant angle-of-attack and flap-deflected cases. For the lift dis- 
tributions due to twist and aileron deflection the presence of the 
fuselage can probably be ignored in most cases. 


Nacelies also tend to affect the lift distribution primarily in 
their own vicinity, but these effects may be significant even for the 
Lift distributions due to twist and aileron deflection. Tip tanks tend 
to increase the lift over much of the outer part of the wing to a large 
extent, particularly in the case of lift distributions due to twist and 
aileron deflection. Except for wings with very high or very low aspect 
ratio, these effects can hardly be underestimated and must be taken 
into account in designing the wing. 


High angles of attack.- Potential flow breaks down at high angles 
of attack, and the higher the Mach number the lower the angle of attack 
at which linearized potential-flow theories such as the one employed in 
this paper fail to predict the Lift distributions accurately. However, 
the critical design loads often occur at high angles of attack. The 
only suggestion that can be made is that once the rigid-wing lift dis- 
tributions at high angles of attack are estimated on the basis of tests 
or experience, the changes in these distributions due to aeroelastic 
effects can be estimated by means of the results calculated in this 
paper. This procedure cannot be justified theoretically because, 
although the nature of the mutual induction effects between various 
parts of the wing after the flow has separated is still substantially 
the same as before, the lifts caused by these induction effects are not 
those predicted by potential-flow theory. However, the changes due to 
aeroelastic action are small unless the speed is near the flutter or 
divergence speed, so that certain inaccuracies can usually be tolerated 
in estimating them. 


longitudinal location of the center of pressure.- The Weissinger 
theory yields no information regarding the location of the chordwise 


center of pressure; however, the assumption that at each spanwise sta- 
tion the section center of pressure of the angle-of-attack loading on 
the two-dimensional airfoil section is unchanged in three-dimensional 
flow has been found by lifting-surface calculations (ref. 6, for example) 
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to be largely justified for swept and unswept wings of moderate and high 
aspect ratio (except near the root and tip). If this assumption is used, 
the longitudinal location of the wing center of pressure may be estimated. 
For low-aspect-ratio wings, the chordwise location of the center of pres- 
gure cannot be determined simply, and lifting-surfece methods must be 
used. For flap and aileron deflections, accurate theoretical methods for 
calculating the longitudinal location of the center of pressure are not 
available, but the approximate methods suggested in reference 7 may be 
applied to obtain qualitative information. 


Effective angle of attack for flap deflection.- In order to deter- 
mine the loading due to flap deflection for wings of high and medium 


aspect ratio (for exemple, А >), the effective angle of attack αρ 


for the flap (or aileron) deflection may be approximated satisfactorily 
by the values obtained from two-dimensional thin-airfoil theory. Fig- 
ure 21 gives a plot of the effective angle of attack ας against flap- 


chord ratio ср/с. For very low aspect ratios (approaching O and cer- 
tainly less than 1/2) values of о, close to 1 are indicated by line- 
arized potential-flow theory, even for relatively small values of сг/с. 


For aspect ratios from about 1/2 to 4, lifting-surface methods must be 
used to obtain potential-flow solutions for the Lift distributions due 
to partial-chord control deflections. 


Calculation of. the roll due to sideslip με The loading for the 


case of full-span ailerons Pail = 1 is the same as the loading on a 


wing with dihedral in yaw or sideslip, because in this case the loading 
on the wing is that due to an angle of attack equal to the product of 
the sideslip angle and dihedral angle on one wing, and the negative of 
that angle of attack on the other wing. The value of Ст [95 for the 


case of full-span aileron deflection is therefore equivalent to the 
parameter Clg 


Some other stability derivatives can be deduced similarly from the 
results presented in this paper. | - 


Relative Merits of the Various Types of 


Aerodynamic Influence Coefficients 


In this paper four types of aerodynamic influence coefficients have 
been discussed: 


(1) The influence-coefficient matrices presented in table V which 
are obtained by solving ος, 8 integral equation by numerical 
methods 
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(2) Influence-coefficient matrices obtained from those of table V 
by multiplying them by interpolating matrices 


(5) Influence coefficients based on Green's function 


(4) Influence coefficients based on flap- and aileron-type lift 
distributions which express the lift distribution in terms of the span- 
wise derivative of the angle-of-attack distribution rather than the 
distribution itself 


The influence coefficients given in table V apply to the stations 
y* = 0.9808, 0.9239, 0.8515, 0.7071, 0.5556, 0.3827, 0.1951, and O. If 
these stations can be used in the calculations in which the influence 
coefficients are to be used, the coefficients given in table V are by 
far the simplest to use because they require no further calculations. 
If these stations cannot conveniently be used, one of the other three 
types of influence coefficients has to be calculated. 


The second type of influence coefficients is based on the first 
and requires a premultiplication of their matrix by a lift interpolating 
matrix and the postmultiplication of that matrix by an angle-of-attack 
interpolating matrix. These interpolating matrices serve to relate the 
lift and angle of attack at the stations of interest to those at the 
stations specified in the preceding paragraph. The interpolating mat- 
rices can be constructed in several ways; parabolic (or possibly cubic) 
interpolation is probably the most satisfactory choice for the angles 
of attack, and for the lift distributions either this type of interpo- 
lation (with a modification at the wing tip) or trigonometric interpo- 
lation should be satisfactory. The interpolating factors do not lend 
themselves readily to automatic computation, but the amount of effort 
involved is relatively small. The two matrix multiplications can then 
be performed readily on automatic computation machines. 


The influente coefficients based on Green's functions are similar 
in concept to the commonly used structural influence coefficients. The 
values of the influence functions [*g'(y*,y,*) and I"&'(y*,yg*) are 


the only aerodynamic influence coefficients discussed herein which indi- 
vidually have physical significance; the first two types of influence 
coefficients have only a collective physical significance in that they 
yield the values of the lift when matrix-multiplied by the sngle-of- 
attack values. However, this individual sigmificance of the coefficients 
based on Green's functions is lost once these coefficients are manipulated 
in the manner indicated in equations (9) and (12) to obtain aerodynamic 
influence coefficients useful in further computations, and the resulting 
influence coefficients have neither more nor less significance than the 
others. The computation of these coefficients requires a relatively 
large expenditure of effort - four matrix multiplications and three mat— 


rix additions, as well as the computation of many values of Г*в’ p οὐ 
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D (by substitution of given values of 8 and 6% in eq. (C3)) 
and a few other minor steps. Despite their conceptual attractiveness, 
these coefficients are therefore practically at a disadvantage compared 
to the other much more simply computed influence coefficients. 


The influence coefficients based on flap- and aileron-type lift 
distributions and on the spanwise derivatives of the angle-of-attack 
distributions are probably the simplest to compute (with the exception 
of those presented in table V); they require only the reading of the 
v&lues of the lift distributions from the figures of this paper &t the 
stations of interest and for the flap and aileron spans corresponding 
to the stations of interest, as well as the multiplication of the mat- 
rix of these coefficients by a diagonal matrix. (For a symmetric dis- 
tribution, a matrix subtraction is also called for.) As previously men- 
tioned, the fact that these coefficients express the lift distribution 
in terms of the derivative of the angle-of-&ttack distribution need not 
be a disadvantage and may actually be an advantuge. The decision as to 
whether to use the second or the fourth type of influence coefficients 
(once the decision has been made that the stations implied in the first 
' type are unsuitable) then becomes largely a matter of individual prefer- 
ence, guided by decisions in amy given case as to the relative conven- 
lence of calculating interpolating factors or reading values from the 
figures in this paper, and of using angle-of-attack distributions or 
their derivatives. 


CONCLUDING REMARKS 


Spanwise lift distributions have been calculated for nineteen 
unswept wings with various aspect ratios and taper ratios and with a 
variety of angle-of-attack or twist distributions, including aileron 
and flap deflections, by means of Weissinger's method with eight control 
points on the semispan. Also calculated by this method were aerodynamic 
influence coefficients which pertain to a certain definite set of sta- 
tions on the span. Three methods for calculating aerodynamic influence 
coefficients for arbitrary stations have been outlined and their rela- 
tive merits discussed. 


The information presented herein can be used in the analysis of 
untwisted wings or wings with known twist-distributions, as well as 
in aeroelastic calculations involving initially unknown twist 
distributions. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., May 5, 1953. 
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APPENDIX A 


MATRIX FORMULATION OF THE WEISSINGER METHOD 
Calculation of the Lift Distribution 


and Influence Coefficients 


From two-dimensional thin-airfoil theory it can be shown that, if 
all the vorticity of a plane or parabolically cambered airfoil section 
is concentrated &t the quarter-chord line, the downwash angle induced 
at the three-quarter-chord line is equal to the geometric angle of 
attack at the three-quarter-chord line. This circumstance leads to the 
Weissinger L-method in which the lifting vortex is concentrated at the 
quarter-chord line and the boundary conditions are satisfied at the 
three-quarter-chord line. The Weissinger equation (ref. 1) can be 
written as 


1 ж 1 x 
ет 3 ar (η) а 1 F(n* аг (n*) ana Al 
a(y*) EL e Cape tirs | (n*,y*) M 7 (A1) 
where 
F(q*,y*) = (пж > о) 
y* - ην 


(m* < 0) 
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and α is the angle of attack or, more specifically, the streamwise 
slope of the mean-camber surface at the three-quarter-chord line. 


Introduction of the trigonometric variables Ө = сов-1 yx and 
$ € cos 1% into equation (Al) yields 


X aux x x 
ale) = = j ar“) а o a F(3,6) аг (9) as (A2) 
hx Јо 48 cos 8 - cos @ 40 dà 


The solution of equation (A2) is effected in reference 1 by the use 
of trigonometric interpolation and integration formulas. An alternate 
solution based on matrix techniques is presented herein which leads to 
the identical results somewhat more simply and suggests a fairly simple 
setup for routine calculations. 


The function Г” is approximated by a finite sine series, as in 
reference 1, so that for any value Э 


r*a = a, sin ng Ἢ | (A5) 
n 


or 


t - [sin оба | o.) | | (A) 


In the calculations of this paper, the values of mn in equation (Ah) 
have been chosen as п = 1, 3, 5, . . . 15 for the symmetrical loadings 
and n=2, 4, 6, . . . Ш for the antisymmetrical loadings; values of ` 


S were chosen at 


шл 


ы ә e. 
n 36 


with m-1,2,5,... 8 for the symmetrical loadings and m = 1, 2, 
52... Т for the antisymmetrical loadings. The use of equal increments 
in ὃ is essential in the method of reference 1, but any values of 3$ 
between O and x/2 could have been chosen in the matrix analysis used 

in the present paper. This possibility of using arbitrary stations in 
this matrix version of Weissinger's method is an important advantage in 
that it allows the direct calculation of influence coefficients for 
arbitrary stations. Such a calculation is complicated, however, by the 
fact that some of the matrices which must be inverted tend to be 111- 
behaved when nonequ&l increments are chosen for 8. ` 
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ar* (a) 


The value of required in equation (A2) can be obtained 


from equation (A5) as 


ar* 
Ew > na, cos ni, (A5) 


so that the first term of equation (A2) becomes 


na, cos nů, dd 
x 2 m 1 sin nêm 
А — ш — > па, ——A 


сов 8 - cos 6, + ^n sin Om 


lx 


Or, in matrix notation, 


The values of a, can be expressed in terms of I", (see eq. (Ah)) ав 


{anb = [sin а |” . А" (A6) 


so that 
x ar* 
ΝΑ ZU cos 8 ET τ) ` ET i 
where 
Be] На аве нарт am) 


The matrices = l, | sin пб |, | sin ое | ^ |, and [Ва | are 


given in table VII(a) for the symmetrical distributions and table VII(b) 
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for the antisymmetrical distributions. As a result of the orthogonality 
of the sine function, the inverse of the [sin пр | matrix is the same 


as one-fourth its transpose except for the first and lest rows. (See, 
for instance, ref. 8.) 


The second term of equation (A2) can be integrated numerically by 
* 
approximating either F or rE by a cosine series. Both approxi- 


mations will yield identical results, and the latter alternative is 
followed here. Thus let 


Е(3,0) ——— аг” το = 2 Ῥηίθ]οος ng 


so that 
π ж 
i | ғ(ә,ө) Г (8) аз = ро(9) (А8) 
Ox Јо d$ 8 


Now, in matrix notation, 


{+ T = [cos nə] ορ (6)* 


so that 


bg(8) = | [eos nôl = 517 лар 
š; |209,9) |; | 


-1 
where | [eos n$| Е is the zeroeth row of the inverse of the matrix 


| cos ne | and. | Τα | consists of the same elements written in the form 
of a diagonal matrix. 
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* 


Inasmuch as the integral of the antisymmetrical component of F = 
is zero, F a ean be written for symmetrical distributions as 
ж x 
^ ag 2 49 
and, for antisymmetrical distributions, as 
А ar, FR + ΕΙ, аг“ 
8 аз 2 49 
where 
Fa(In*y*) = r(5,y*) (n* > o) (A9) 
Fr(|n*|,y*) = F(n*,y*) (1% <o) (A10) 


can be expressed in matrix 


From equations (A5) and (A6) a 


notation as 


T8 e Eos Πο ай} 


so that equation (A8) can be written for symmetrical distributions as 


E IN Fa (8, να dò ΕΙ A [cos nën | |n] [sin 2181] ES Ta 
TNEGA an) 


where 


Pa = ыы Eo =] sn]? 
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Similarly, for antisymmetrical distributions 


E ο ва] «ehh ip waa) 
where “7 ү 
BJ «dsl tf e? 


The matrices [D], [ту |], апа. | сов nox | п] are given in table VII (a) 


for symmetrical distributions and table VII (b) for НЫ 
distributions. 


Equation (Al) (or its equivalent, eq. (A2)) can now be expressed 
completely in matrix form. For symmetrical distributions, the equa- 


[5η ή - 63 
E t - 4%) (A15) 


and for antisymmetrical distributions, 


[δε] - eF] [Pa] ra P {aa} 


or | E m 


ALAS {a,b (АЈА) 


It should be noted that the [в] end [D] matrices are invariant 


or 


with plan form and that only the [F] matrices need be computed sepa- 
rately for each plan form; ali the matrices are independent of the 
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angle~of-attack conditions. A computing form for the elements of the 
F matrices is given in table VIII; this computing form includes pro- 


vision for calculating the load on swept wings. Sample 2 [Fa] and ДЕЙ 


matrices аге shown in table IX. 


Equations (A15) and (А1%) can be expressed as 
t5 = [e] a) (19) 


so that the elements of the matrices [a] 7l constitute, in effect, sets 


of aerodynamic influence coefficients. The influence coefficients pre- 
sented in table V for the plan forms treated in this paper are defined 
as 


= -i-[a, == 
Bl = 2-4 
and 
Е -1 
[ος = zz 
so that 


4] = сов} 


and (A16) 
ы) т ACIES: 


The division by m and by ia has been performed both to facilitate 


interpolation of the coefficients for unswept wings with plan forms other 
than those considered in this paper and for convenience in aeroelastic 
calculations. Inasmuch as the lift distribution is much less sensitive 
to Mach number than is the over-all magnitude of the lift, an influence- 


coefficient matrix [3 or - [aa] chosen for the average of the sub- 
sonic Mach number range of interest (that is, for the effective aspect 
ratio ΑΝΙ - M? corresponding to that average Mach number) will serve 
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for the entire range, provided only that for each Mach number the appro- 
priate values of Cr and ση. are used in equations (А16). (See 


ref. 9, for instance, for simple methods of estimating Mach number 
effects on Cr and C, ) 
d. 


Calculation of the Aerodynamic Parameters Associated 
With the Lift Distributions 
Ihe values of the lift, induced-drag, bending-moment, and rolling- 
moment coefficients can be obtained conveniently by the use of the inte- 
grating matrices derived in this section. 


An integrating matrix for the lift coefficient associated with 
symmetrical loadings can be obtained as follows: ` 


The lift coefficient can be written as 
1 
2/0 


If, ав before, у* = сов Ө and P. = > 8, Sin пб, with n = l, 5, 
n 


Б, < « « №, then 


A x/2 | Е 
ст, = an > an sin nO sin Ө 49 
n 
or . _ 


where 
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and | [sin a8] | is the first row of the matrix ‘sin пен“ given 
l 


in table VII(a); the matrix | Тоу | is given in table X. 


Similarly, the integration to obtain the bending-moment coefficient 
for symmetrical loadings Съм can be performed as follows: 


The bending-moment coefficient can be written as 


1 


- A(2 τα. ode e эй od S 
= 8 t 595 "δι 85 T 15 9 7" тт * 
—2 a -a +22 ) 
117 ++ 165 12 ош 15 
Or 


om = роза 4779) 


where the mstrix 


1|2 2 2 2 2 2 2 2 -1 
=== = 2-5 m -Ξ- --- --Ξ-  — ып no 
Б C] ΠΕ 5 21 15 77 117 165 | [sin ne] 
is given in table X. 


For antisymmetrical loadings, the half-wing lift coefficient “hy ур 


is obtained in the following manner: 


The equation for the half-wing lift coefficient is 


A ш * 
туа EJ, Гау” 
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and, if p = > a, Sin n0, with n = 2, ‚б, . . . Ш as before, 
> | 


х/2 | 
ay 
= = a, sin nð sin Ө аө 
“Ly Ja о 2. n š 


= А2 а - & Dared 2. . -6- == 
A (5 вв = 35 6 = 25 Sp * LOL E. ED 


ог 


“Ly уә > la, Қ TD 


where the matrix 


8 


is given in table X. 


lil 
R 
[| 
| o 
Ivi 
I 
| 
Ul 
ON 
--3 
Го | 
5. 
B 
5 
[|] 
ЕЗ 


Similarly, the rolling-moment coefficient can be obtained as follows: 


The equation for the rolling-moment coefficient is 


A ы ж 
а-%/ ΓῈ уж дуж 


x/2 - 
-EJ > a, Sin пе cos Ө sin Ө ἄθ 
О п 


ll 
Bla 


or 


ο 
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where the matrix 


is given in table X. 


The induced-drag coefficient Ор, can be written ав 


where 
41-42% 
1 аг s(7*) x 
πο = 2 sn mmt. 
Вх 1 44% y* - ņ* 
or 


iy ERIS 


For symmetrical loadings the τοι] matrix can be used to integrate 


the values of o DT. go that 


°p, = alto, 4а) 


An integrating matrix to evaluate Cp, for antisymmetrical distributions 
can be set up similarly. However, in this paper, values of Ср. were 


calculated only for additional lift distributions. 
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APPENDIX B 


CALCULATION OF LIPT DISTRIBUTIONS FOR DISCONTINUOUS 


ANGLE-OF-ATTACK CONDITIONS 


The method of solving equation (Al) outlined in appendix A relies 
heavily on numerical integration, as does the method of reference 1. 
Discontinuous angle-of-attack distributions therefore cannot be analyzed 
as accurately as continuous ones can, because discontinuous angle-of- 
attack distributions are known (on the basis of knowledge of the lift dis- 
tributions of wings with very low and very high aspect ratios presented in 

x 

ref. 3) to give rise to logarithmic singularities in the function т. 

T 
which occurs in the integrands of both integrals in equation (Al). Nor 
can a discontinuous angle-of-attack distribution be described adequately 
by a small number of points on the semispan; for instance, with stations 
located as they are for the calculations described in this paper any 
inboard flap terminating at a value of y* greater than 0.3827 but less 
than 0.5556 would, for a unit effective angle-of-attack distribution, be 
characterized by the angle-of-attack distribution 1, 1, 1, O, O, O, O, O 
regardless of the exact location of the end of the flap. 


These difficulties can be overcome by using the results obtained by 
solving equation (Al) for the case of wings of vanishingly small aspect 
ratio (see ref. 3); in which case the second integral vanishes. This 
technique is similar to the one used by Multhopp (ref. 10) in connection 
with the Prandtl lifting-line equation to handle discontinuous angle-of- 
attack distributions. The lift distribution Г* is considered to consist 


of & "discontinuous" part г” which is the solution to equation (А1) if 


the second term is neglected and of & correcting "continuous" part γα, 
so that 


Г“ = P*a + p. (B1) 


where r^, is defined implicitly by 


т * 
ar 
o 48 cos 3 - cos 8 
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and г“ is the correction that must be added to Γη to obtain 8 


function Г* which satisfies equ&tion (Al) for the given discontinuous 


angle-of-attack distribution. The solution of equation (82) for r" 


corresponding to the more common discontinuous angle-of-attack distri- 
butions is given in reference 3; specifically, for inboard flaps termi- 


nating at y* = Yo* = совт 96 


. 3 + Bo 
Ы l sin —5— 
Г p, 3590) - τ] - 205)sin 3 - (cos 3 - cos θο)198ς — Te = e| = 
sin 2 
2 
ὃ + 8o 
cos 2 
(сов 9 + сов 05)1oge (вз) 
ὃ - Ө 
cos 
2 
and for outboard ailerons with inner ends at y* = y-* = cos-i θο 
ὃ + 9g 
x lh sin —— 
Г Dai (9:90) = = (cos 9 - сов θο) loge — 1$ - 85) = 
Sin — 
2 
ð + 6 
сов —— 
(cos 8 + cos 65)1og, ———— — (ві) 
8 - Og 
сов ------ 
2 


If equations (Bl) and (B2) are substituted into equation (A2), the 
result. is 


Ы ж 
ал О 99 сов 9 - cos Ө Әл“ о ( 2 ) 39 ( ? o) 


(85) 


4O NACA TN 3014 


where the function R(8,0,) is defined by 


R(0,0,) = L pg) 98и - A a 


Or, specifically, for flaps 


ZODEESAM Е. (9,8) 


and for ailerons 


R.(9,9o) = 2 |" ғ.(9,0) > ai (9) as (B6b) 


Comparison of equation (B5) with equation (A2) indicates that equa- 
tion (B5) may be considered to be the Weissinger equation (eq. (Al) or 
(A2)) for the lift distribution Г*с оп the given wing (the plan form 


of which determines the function Р(9,6)) corresponding to ап angle-of- 
attack distribution R(0,0,). Inasmuch as к(ө,6.) ів а continuous 


function, as is demonstrated presently, equation (B5) сап be solved in 
the manner used for equation (A2). If (0,02) is being evaluated at 


the stations considered in this paper ( = E 2 ao. 3 E e. . 5), апа 
16 16 ' 16 e 


if the eight values of R are listed in & column in the order of 
increasing Ө, then premultiplication of this column by the matrix [a] 
given in this report and by the appropriate value of s yields the 
desired function P*a for the given discontinuous angie-of-attack 
distribution. 


As indicated in equations (86), the function R(0,05,) depends on 


the plan form, which determines Ұ(9,0), and on the position of the dis- 
continuity in the angle-of-attack distribution, which determines г” 


ar * (9,8) | 
3 ` = 


and, hence, 


6A 
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For flaps and ailerons, 


9 + 8, + 96 | 
Зг"р, " Een cos — 
= — (π - 28,)cos 8 + sin ð| log, + log, 
99 к |8 - 8o | 9 - 96 
sin сов 
2 2 | 
(B7a) 
and 
9 + 6 ὃ + Ө 
EM ji сов —— —— sin —— 
————— = —/2 sin 9. + sin 39110 - lo B 
T = 5 Ec s 8 82 ΠΕΝ (BTb) 
сов ---------- Sin-————— 
2 2 


ar* 


во that the functions may be seen to have logarithmic singulari- 


ties. The evaluation of R(8, 8, ) from equations (B6) by numerical 


methods is therefore not a trivial problem. A logarithmic singularity 
is integrable, however, and F(9,0) is always continuous in 9 and Ө 
so that R(9,05) must always be continuous. The integration can thus 


be effected readily by expanding F($,0) in а finite Fourier series as 
follows: 


Let 


Р.(9,0) = > Pa (9)сов а- (а ео) 
n 


(88) 


Εα(8»θ) -=> Pg (0)cos nd (n=O, 2, Ше. ο αθώο) 


n 


Substitution of these expressions for F and those of equations (BT) 
x 


ar 
for T into equations (B6) yields 


Rs(9,9o) = 25 У в, (9) [( - 98ο)τα + Jn + Καὶ 


(в9) 


Βα(θ»θο) = 2 Ps, (7a n Kn) 
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where 


л 
In = | сов 9 cos nd 4% 
ο ; 


9 + 8 
π sin 5 
Jn = | sin 9 cos nd loge — —l18 - 6.T а 
0 ҒАЛЫ (810) 
2 
3 + 6, 


π σος — _ 
Kn = | sin ð cos n8 log, —— ————— 49 
O 9 - 8o 
сов -------- 
2 


These integrals сап all be evaluated explicitly and аге 


Io = 9 
= £ 
Ту = 5 (в11) 
In = 0 (n= 2, 5,...) 
Jo = x sin 9g 
J, = ® sin 29 
La we о (B12) 
д|вїп(п + 1)θο sin(n - 1}θο 
J, = > n = 2 "Er 
n i ne 1 n= l ( ‚ 2» ) 
and 
Ko = -x sin 8, 
π 
Κι = — sin 20 
ΕΣ 5 (із) 
nl 4 8in(n + 1);  sin(n - 1), 
K, = (-1) Ὁ di a SN ERU RS CNN (n = 2, 3, xx x) 


n + 1 n - 1 
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With these values for Ip, Jp, and Κι, equations (B9) сап be simplified 
to 


RC (0,86) = = У ;ερ(θ»θο)Βα(ϑο) (n= 1, 5, 5, + + + 19) 
n 
(він) 
Βα(θ»θο) = => Ps (8285)En(8o) (n = 0, 2, by... 16) 
n 
where 
Ho = 0 
_ x - 20o + sin 265 
E 2 (B15) 
sin(n + 1)9 sin(n - 1)e 
H = о _ о (n & δ, 5, + 4) 


m + 1 n - 1 


The function R(0, ео) can also be expressed in matrix form by writing 
equations (B8) in matrix form as 


|F (3,6) | ESOS] [cos n8] 


lil 


so that 


үшү 
| (9,0 )| [eos " (B16) 


[sot 


! 
where | cos n$ | is the transpose of the matrix | cos nð]. The following 


expressions are then obtained by combining equation (B16) with the matrix 
equivalent of equations (Bl): 


в (0,85) = dk raso) | [6ος 2317] ^ ве), 


1 TED Г ü (B17) 
Ra (0,90) m [re (а,в) [Гов nd] | qFn(9o) | 
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In essence the procedure used in this appendix for calculating lift 
distributions for discontinuous angle-of-attack distributions can thus 
be seen to consist in avoiding the necessity of integrating numerically 
an initially unknown singular function, as required in equation (Al) or 
(A2), by performing, in effect, the integration of the singular part of 
the function analytically (by solving eq. (B2), as in ref. 3) and 
treating only the continuous part of the function numerically (by 
solving eq. (B5)). As part of this procedure, a singular function has 
to be integrated numerically in order to evaluate Ἀ(θ;θϱ); however, 


this function is known initially because ít is the product of a known 
singular but integrable function and a known continuous function 


ж 
T' 
| — and Р(9,0), respectively] so that, by expanding the continuous 


function in a Fourier series, the numerical integration can be effected 
without difficulty. `. 


-1 
1 
The values of | [eos ná] | ENOS for flap-span ratios of 0.1, 


0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9 and for aileron-span ratios 
of 0.1, 0.2, 0.5, O.4, 0.5, 0.6, Ο.Τ, 0.8, 0.9, and 1.0 are given in 
table XI. 
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APPENDIX С 
CALCULATION OF GREEN'S FUNCTIONS 


As pointed out in the body of this paper (see eqs. (3) and (№)) 
the desired Green's functions are 


ar*e( yy," ) _ СИТИ 
ду” дуо“ 


The functions га and PES can be obtained in the manner indicated 


in appendix B. However, in order to calculate the desired derivatives 
of these functions, & numerical differentiation would have to be per- 
formed with respect to yg*- Sueh & numerical differentiation is inher- 


ently inaccurate inasmuch as po and pog have singularities. The 
desired Green's functions are therefore best calculated without using 
the calculated values of Ге and Pg. and by using, instead, & 
modification of the method of appendix B. 


Let 
3r" (уу) 
oy,* 


г”; (уу) = 
(C1) 
X 
or ai1(s*s¥o*) 
Ον. 


1 
Then p and P. are lift distributions corresponding to impulse- 
type angle-of-attack distributions of the following forms: 


D (y =- 


O. 


Symmetrical loading Amtisymmetrical loading 
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Hence, they must satisfy equation (Al) or its equivalent, equation (A2), 
for these angle-of-attack distributions. 


Again, as in appendix B, the lift distributions can be considered 
to consist of a discontinuous part, which satisfies equation (Al) for 
the given angle-of-attack distribution if the second integral on the 
right side is disregarded, and of a correction part; that is, 


x ! x ! p* | 
8 p S с 
(C2) 


t 1 
The functions I, р and Га р must thus satisfy equation (B2) for 


the given angle-of-attack distributions. By virtue of the linearity of 
equation (B2) and by virtue of the definitions of the functions in terms 


| _ t? 
of Г" апа ЕРТ respectively, prs Б апа Ps can be obtained 


by differentiating with respect to yg" the solutions of equation (B2) 


given in equations (B5) and (Bh) for the flap and aileron angle-of- 
attack conditions. Thus 


' sin 9 + sin Ө 
Ы | sin ð - sin θο | 
(C3 ) 
! sin(9$9 + Ө 
віп |9 = θο | 
Similarly, the functions Is с and I шив satisfy equation (B5), 


where now В (е, 8 5) is defined for the symmetrical and ΕΞ 


loadings, respectively, by 


π ους (9,0 
Rs (9,95) = E Fg (8,0) -- as 


and | (Ch ) 


ara (828) 


X 
R (6,05) = & J. F, (3,8) — ë 
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The evaluation of these integrals can be effected in the manner 
employed for equations (B6) so that 


Βρ(θ»θο) = ж2 Pen(®)E'n(80) (n= 1, 3, 5, + + + 15) 
(C5) 
Ва(9,6о) = => Ps (0)E'n(80) (n= 0, 2, à, . . . 16) 
where 
H π(θο) = 2 sin πθο (C6) 


for all values of n, and where Pan and Р. are the same values as 
n 


those used in appendix B. Thus, in matrix form, for given values of 0 
and 6 
о? 


Βε(θ,θο) = БЖС (Гео a {sin neob 


(CT) 


Βα(θ»θο) = Lip, (8,0)] | [eos ΙΝ = {sin пёс} 


r] -1 
The values of | [eos 28] | {sin ne, b are given in table VI for 


values of Уы” ranging from O to 0.9 for symmetrical distributions and 
0.1 to 0.9 for antisymmetrical distributions. 

The desired Green's functions can thus be calculated in the fol- 
lowing way. For a given value of Өс, the values of Rg(9,09) and 
Βρ(θ;60) are calculated for eight equal increments of Ө between О 


and х/2 from equation (CT). These values are then written ав 
columns and premultiplied by the matrix of influence coefficients 
tabulated in this paper for the given plan form in order to obtain 


в (9,00) Ра (8,89) 


στο Cla 
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for the given value of 9-. То these values are added the values of 


f 
P^. (9,9) -— г“ (8,80) 
Cr, Cza 


obtained by dividing the values of Pa 5 and Pa calculated from 
equation (C4) by On, and e respectively, for this value of Ө, 
and the given values of Ө. This procedure yields ` 


r*a (0,80) а Га (6,99) 
Cr Ста 


( The division by Ст, and by uM is performed to facilitate the 
further calculations required to obtain the desired influence coeffi- 
cients, as explained in the body of this paper.) This calculation is 


repeated for all the values of 98, for which the Green's functions are 
desired. 


TA NACA TN 3014 49 


- REFERENCES 


т l. Weissinger, J.: The Lift Distribution of Swept-Back Wings. МАСА 
TM 1120, 1947. 


2. Diederich, Franklin W.: Calculation of the Aerodynamic Loading of 
Swept and Unswept Flexible Wings of Arbitrary Stiffness. NACA 
Rep. 1000, 1950. 


5. Diederich, Franklin W., and Zlotnick, Martin: Theoretical Spanwise 
Lift Distributions of Low-Aspect-Ratio Wings at Speeds Below and 
Above the Speed of Sound. NACA TN 1973, 1949. 


4. Diederich, Franklin W.: Calculation of the Lateral Control of Swept 
and Unswept Flexible Wings of Arbitrary Stiffness. МАСА Rep. 1021, 
1951. 


5. DeYoung, John, and Harper, Charles W.: Theoretical Symmetric Span 
loading at Subsonic Speeds for Wings Having Arbitrary Pian Form. 
NACA Rep. 921, 1948. 


6. Falkner ; У. M.: The Calculation of Aerodynamic Loading on Surfaces 
of Any Shape. В. & M. No. 1910, British A.R.C., 1943. 


T. Diederich, Franklin W.: A Simple Approximate Method for Calculating 
Spanwise Lift Distributions and Aerodynamic Influence Coefficients 
at Subsonic Speeds. NACA TN 2751, 1952. 


8. Benscoter, Stanley U.: Matrix Development of Multhopp's Equations 
for Spanwise Air-~Load Distribution. Jour. Aero. Sci., vol. 15, 


" 9. Diederich, Franklin W.: A Plan-Form Parameter for Correlating Certain 
Aerodynamic Characteristics of Swept Wings. МАСА TN 2555, 1951. 


10. Multhopp, H.: ‘The Calculation of the Lift Distribution of Aerofoils. 
R.I.P. Translation No. 2592, British M.A.P. 


20 МАСА TN 301^ 


TABLE I 


INDEX TO SPANWISE-LIFT-DISTRIBUTION FIGURES 


[А11 wings аге unswept | 
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TABLE IT 


FORCE AND MOMENT COEFFICIENTS FOR ANGLE-OF-ATTACK LOADINGS 


ЕТТЕ 


| Additional loading — loading Damping-in-roll loading 


ЕЛСЛЕЛЕЛЕСЕК OE 


2. 
5% 
3: 
5. 
5s 
1. 
Ц. 
4, 
1. 
h. 
5. 
5. 
5. 
5. 


TABLE III 
LIFT AND BEADING-MOMERT COEFFICIENTS FOR UNIT EFFECTIVE FLAP DEFLECTION 


ЕР o2 I. oA Ἔ κο) Жо. 2.0.9 
Plan form} Р 


τ πι. m | % 


«33186 1.217h6| . 
| .25hh5|1.01760| . 
53222 1.2004 . 
‚22729 1.17816| . 


—— 


„#156 |1.92581 | . 
| „50260 1.955061] . 
„5601611..57482| .50225 1.95017 |. 
. 35166 [1.51335 9080 1.85973 Е 
.21417|1.09860| .55891 |1.95065 | .47517/1. 76911) . 


eG 


πτος NIL VOVN 


TABLE IV 
LUFT AND ROLLING-MOMENT COEFFICIENTS FOR UNIT EFFECPEIVE AEROS DEFLSCTION 


πτος NL VOVN 


Dec De [s wes Des eo Doe De 
Plan form 
eS Sos eae eS SSS 


~6760|1.26012| . 
-63168/1.41156] . 
.69120[1.55088| . 


eG 


0.0530 
0.0224 
0.0121 
0.0095 
0.0075 
0.0068 
0.0064 
0.0065 


looooooool 


0.0590 
0.0274 
0.0154 
0.0116 
0.0091 
0.0080 
0.0073 
0.0072 


looooooool 


0.0591 
0.0277 
0.0157 
0.0119 
0.0095 
0.0081 
0.0073 
0.0073 


[оооооооо! 


0.0283 
0.0162 
0.0122 
0.0094 
0.0081 
0.0072 
0.0070 


[oooooooo! 


0.0599. 


TABLE V 
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AERODYNAMIC-INFLUENCE-COEFFICIENT MATRICES 


оооооооо 
о 
VN 
= 
H 


ОООООСОСОСОО. 
о 
= 
2] 
NO 


оооооооо 
NO 


оооооооо 
ON 
Со 


оооооооо ' 
2 
E 


(a) Symmetric loadings Е 


Plan form 311. 


.0218 
0644 
«1525 
. 5275 
‚1951. 
.1426 
.1297 
«1522 


оооооооо 


Plan form 512 


оооооооо 
VN 
№ 
Ον 
NO 


Plan form 313 


ovo 
. 0893 
«1709 
5292 
1971 
1h26 
.1276 
1206 


оооооооо 


Plan form 514- 


0.0459 
0.0925 
0.1752 
ο. 3344 
0.1991 
0.1419 
0.1248 
0.1169 


0.0202 0.0182 0184 0.0086 
0.0587 0.0517 0.0508 0.0245 
0.1144 0.1041 0.0956 0.0475 
0.2221 0.1722 0.1614 0.0768 
0.4276 0.2979 0.2452 0.119} 
0.2746 0.5320 0.5910 0.1720 
0.2170 0.3737 0.6582 0.2788 
0.2114 0.5258 0.5497 0.4974 
0.0556 0.0558 0.0518 0.0159 
0.0782 0.0704 0.0680 0.0550 
0.1296 0.1189 0.1095 0.05hÀ 
0.2269 0.1789 0.1676 0.0801. 
0.1270 0.2929 0.2599 0.1178 
0.2676 0.5115 0.5759 0.1655 
0.2106 0.3581 0.6274 0.2656 
0.2047 0.3118 0.5253 0.4733 
0.0502 0.0565 0.0343 0.0171 
0.0826 0.0747 0.0721 0.0550 
0.1342 0.1256 0.1140 0.0566 
0.2300 0.1821 0.1706 0.0815 
O.4174 0.2955 0.2400 0.1178 
0.2662 0.5083 0.3725 0.1637 
0.2079 0.3536 0.6209 0.2622 
0.2014 0.5066 0.5161 0.4688 
0.0401 0.0382 0.0561 0.0179 
0.0860 0.0780 0.0752 0.0565 
0.138 0.1274 0.1174 0.0582 
0.2341 0.1855 0.1751 0.0826 
0.4215 0.2948 0.2400 0.1176 
0.2654 0.5083 0.3695 0.1616 
0.2035 0.3480 0.6160 0.2588 
0.1955 0.2986 0.5077 0.4664 
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0.0609 
0.0288 
0.0165 
0.0124 
0.0095 
0.0080 
0.0070 
0.0067 


оны 


0.0272 
0.0094 
0.0042 
0.0051 
0.0025 
0.0020 
0.0018 
0.0019 


Гоооооооо!| 


0.0159 
0.0080 
0.0056 
0.0040 
0.00354 
0.0030 
0.0029 


Гсооооооо! 


0.0367 
0.0165 
0.0086 
0.0060 
0.0042 
0.0055 
0.0050 
0.0029 


0.0565 


IABLE V.- Continued 


AERODYNAMIC-INFLUENCE-CORFFICIENT MATRICES 


0.0565 
0.1519 
0.0808 
0.0509 
0.0400 
0.0222 
0.0288 
0.0270 


.0128 
0716 
.0509 
0164 
.0128 
. 0102 
. 0096 
.0092 


ΟΟΟΟΟΟΟΟ 


0.0509 
0.0864 
0.0414 
0.0232 
0.0173 
0.0136 
0.0124 
0.0118 


0.0525 
0.0887 
0.0457 
0.0248 
0.0182 
0.0110 
0.0124 
0.0117 


0.0470 
.117һі 
.21:75 
. 1399 
.0933 
. 0761. 
.0651 
.0627 


ΟΟΟΟΟΟΟ 


. 0087 
. 0102 
1245 
.0607 
.0565 
.0502 
.0262 
. 0261 


ОООООООО 


0.0219 
0.0591 
«1252 
.0688 
.Ol1T7 
‚0552 
.0285 
.0281 


OOOOOO 


.0O2h2 
. 0631. 
‚1595 
.0720 
O45) 
0340 
. 0284 
.0277 


ΟΟΟΟΟΟΟΟ 


(a) Continued 


Plan form 415 


0452 
0945 
. 178} 
. 5598 
.2011 
„1%15 
„122} 
„1152 


ΟΟΟΟΟΟΟΟ 


0.0410 
0.0878 
0.1408 
0.2374 
0.4261 
0.2655 
0.1999 
0.1900 


Plan form 501 


0.0062 
.O2hl 
.0718 
.1795 
.0961 
0641 
0569 
0534 


ΟΟΟΟΟΟΟ 


. 0061 
. 0209 
.0467 
.lOTT 
‚2555 
. 1385 
.1029 
. 1000 


OOOOOOOO 


Plan form 322 


0.0186 
0.0404 
0.0855 
0.1820 
0.0985 
0.0648 
0.0564 
0.0525 


0.0150 
0.0340 
0.0584 
0.1127 
0.2281 
0.1528 
0.0970 
0.0936 


Plan form 423 


осооооооо 
5 
Fa 


0.0171 
0.0579 
0.0650 
0.1166 
0.2296 
0.1517 
0.0942 
0.0901 


0.0390 
0.0795 
0.1293 
0.1871 
0.2964 
0.5104 
0.3439 
0.2912 


0.0047 
0.0174 
0.0505 
0.0746 
0.1483 
0.2945 
0.1911 
0.1647 


0.0135 
0.0285 
0.0501. 
0.0795 


0.1441 


0.2774 
0.1801 
0.1512 


0.0152 
0.0314 
0.0557 
0.0820 
O. LHi 
0.2746 
0.2.717 
0.11448 


0.0367 0.0182 
0.0764 0.0370 
0.1187 0.0588 
0.1740 0.0828 
0.2597 0.1111 
0.3677 0.1600 
0.6158 0.2564 
0.5005 0.4657 
0.0057 0.0022 
0.0169 0.0080 
0.0555 0.0177 
0.0615 0.0515 
0.1117 0.0549 
0.2011 0.0850 
0.5696 0.1520 
0.3014 0.2910 
0.0122 0.0061 
0.0266 0.0126 
0.0437 0.0216 
0.0709 0.0550 
0.1080 0.0525 
0.1874 0.0788 
0.5427 0.1400 
0.2772 0.2703 
0.0156 0.0067 
0.0291 0.0138 
0.0565 0.0228 
0.0727 0.0557 
0.1072 0.0517 
0.1851 0.0764 
0.5552 0.1557 
0.2682 0.2655 


56 


[ooooo0o000l ΣΕ [Dooooooo! [Doooooool 


ΟΟΟΟΟΟΟΟ 


оооооооо« 


оооооооо 


оооооооо: 


TABLE V.- Continued 
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AERODYNAMIC-INFLUENCE-COEFFICIENT MATRICES 


oooooooo 


оооооооо 


(а) Continued 


Plan form 424 


| 0.0252 


0.0489 
0.0965 
0.1955 
0.1036 
0.0651 
0.0534 
0.0881 


0.0190 
.0415 
.0678 
.1218 
.2551 
„151% 
. 0900 
. 0813 


OOOOOOOQ 


Plan form 225 


0.024} 
0.0513 
0.1005 
0.2000 
0.1065 
0.0652 
0.0517 
0.0455 


0.0200 
0.0435 
0.0706 
0.1257 
0.2409 
0.1318 
0.0868 
0.0795 


Plan form 551 


. 0010 
. 0076 
0515 
.1059 
ο η 
. 0249 
. 0212 
.0191 


оооооооюо 


.0019 
.0061 
.0165 
.0500 
137% 
. 0618 
. 01.26 
„0422 


OOOOOOOO 


Plan form 432 


0.0080 


0.0175 
0.0432 
0.1090 
0.0482 
0.0258 
0.0210 
0.0186 


0.0055 
0.0152 
0.0255 
0.0548 
0.1340 
0.0642 
0.0395 
0.0577 


оооооооо 


OO000 000 


0.0147 0.0072 
0.0510 0.0156 
0.0586 0.0257 
0.0740 0.0340 
0.1061 0.0507 
0.1785 0.0722 
0.5285 0.1506 
0.2562 0.2612 
0.0150 0.0073 
0.0516 0.0147 
0.0892 0.0238 
0.0745 0.0557 
0.1050 0.0497 
0.1756 0.0707 
6.3243 0.1267 
0.2461 0.2586 
0.0017 0.0003 
0.0045 0.0022 
0.0110 0.0055 
0.0246 0.0109 
0.0462 0.0222 
0.1016 0.0598 
0.2192 0.0846 
0.1693 0.1850 
0.0040 0.0020 
0.0089 0.001 
0.0187 0.0073 
0.0261 0.0115 
0.0434- 0.0206 
0.0913 0.0346 
0.1985 0.0748 
6.1495 0.1688 
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0.0263 
0.0108 
0.0048 
0.0029 
0.0017 
0.0013 
0.0010 
0.0010 


0.0275 
0.0118 
0.0055 
0.0055 
0.0020 
0.0014 
0.0010 


0.0010 


0.0287 
0.0124 
0.0059 
0.0036 
0.0021 
0.0014 
0.0010 
0.0009 


TABLE V.- Continued 


AERODYNAMIC-INFLUENCE-COERFFICIENT MATRICES 


0.0211 
0.0602 
0.0258 
0.0120 
0.0077 
0.0051 
0.0045 
0.0059 


0.0251 
0.0644 
0.0289 
0.0138 
0.0086 
0.0054 
0.0045 
0.0037 


0.0244 
0.0677 
0.0310 
0.0150 
0.0092 
0.0055 
0.0041. 


0.0025 


0.0133 
0.0372 
0.0895 
0.0388 
0.0188 
0.0151 
0.0098 
0.0096 


0.0156 
0.0420 
0.0968 
0.0552 
0.0209 
0.0158 
0.0096 
0.0090 


0.0169 
0.0451 
0.1025 
0.0465 
0.0222 
0.0140 
0.0092 
0.0083 


(а) Continued 


Plan form 333 


0.0100 
0.0215 
0.04.87 
0.1140 
0.0505 
0.0263 
0.0206 
0.0179 


0.0069 
0.0159 
0.0271 
0.0585 
0.1359 
0.0655 
0.0374 
0.0350 


Plan form 551 


0.0120 
0.0256 
0.0550 
0.1220 
0.0540 
0.0268 
0.0196 
0.0163 


0.0083 
0.0188 
0.0312 
0.0655 
0.171 
0.0635 
0.0345 
0.0511 


Plan form 555 


0.0091 
0.0205 
0.0338 
0.0672 
0.1468 
0.0637 
0.0521. 
0.0281 


0.0057 
0.0115 
0.0206 
0.0550 
0.0691 
0.1582 
0.0833 
0.0618 


0.0065 
0.0141 
0.0229 
0.0451 
0.0705 
0.1582 
0.0773 
0.0555 


0.0069 
0.0159 
0.0211 
0.0363 
0.0717 
0.1597 
0.0752 
0.0475 


0.0056 
0.0101 
0.0160 
0.0269 
0.0425 
0.0870 
0.1906 
0.1394 


0.0051 
0.0110 
0.0169 
0.0272 
0.0407 
0.0820 
0.1824 
0.1254 


0.0052 
0.0112 
0.0170 
0.0271 
0.0594. 
0.0786 


0,1777 


0.1156 


0.0025 
0.0046 
0.0078 
0.0117 
0.0199 
0.0522 
0.0702 
0.1651. 


ΕΓ. 


0.0025 
0.0048 
0.0081 
0.0115 
0.0187 
0.0290 
0.0640 
0.1565 


ЕЕ 


0.002 
0.0048 
0.0080 
0.0111 
0.0177 
0.0265 
0.0594 
0.1520 


Plan form 341 


0.0075 -0.0008 0.0015 -0.0003 0.0008 -0.0002 0.0007 -0.0001 


0.0008 
0.0002 
0.0002 
0.0001 
0.0001 
0.0001 
0.0001 


0.0250 
0.0044 
0.0012 
0.0010 
0.0005 
0.0006 
0.0005 


0.0040 
0.0422 
0.0109 
0.0035 
0.0028 
0.0019 
0.0021 


0.0025 
0.011) 
0.0631 
0.0196 
0.0077 
0.0064 
0.0052 


0.0012 
0.0047 
0.0207 
0.0845 
0.0508 
0.0149 
0.0148 


0.0014 
0.0056 
0.0089 
0.0519 
0.1062 
0.0471 
0.0516 


0.0009 
0.0029 
0.007 
0.0159 
0.0476 
0.1545 
0.0918 


0.0006 
0.0015 
0.0051 
0.0077 
0.0157 
0.0452 
0.1247 
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. 0195 
‚0059 
.0017 
. 0009 
.0005 
‚ 0004. 
. 0003 
. 0005 


.0208 
.0072 
. 0021 
. 0012 
.0006 
‚ООО 
.0003 
.0005 


.0222 
. 0084 
. OO21 
.0016 
0007 
0005 
0003 
0005 
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TABLE V.- Continued 


NACA TN 3014 


AERODYNAMIC-INFLUENCE-COEFFICIENT MATRICES 


0113 
0387 
.0115 
. 0036 
. 0022 
. 0013 
. 0012 
. 0010 


. 01.1. 
. 0456 


0146 


0049 
. 0028 
.0015 
. 0015 
.0011 


. 0165 
. 04Ł84— 
.0178 
. 0061. 
. 003} 
. 0017 
.0013 
.0010 


oOo00 0000 


O9O000000 


оооооооо, 


. 006 


0161 


‚ОБЕ 
. 0165 
.0056 
.0038 
.0025 
.0028 


.0066 
.0210 
.0602 
.0195 
.0067 
.ooke 
.0026 
. 0027 


. 0088 
«0259 
. 0675 
.0252 
. 0081 
.0056 
. 0027 
. 0025 


(a) Concluded 


Plan form 342 


0.0029 
.0062 
.0198 
.0688 
.0219 
. 0083 
0064 
.0051 


ооооооо 


.0018 


.0078 
. 0216 
.0857 
.0290 
‚0152. 
.0126 


ΟΟΟΟΟΟΟΟ 


Plan form 343 


0.002 
0.0088 


O. 024}-- 


0.0736 
0.0237 
0.0087 
0.0063 
0.0049 


Plan form 34h- 


0.0057. 


0.0118 
0.0298 
0.0805 
0.0262 
0.0091 
0.0060 
0.0045 


0.0025 
0.0057 
0.0097 
0.0275 
0.0855 
0.0286 
0.0120 
0.0115 


0.00317 
0.0071. 


0.0121 
0.0510 
0.0896 
0.0281; 
0.0107 
0.0095 


. 0013 


оооооооо 


оооооооо 


OOOOOOOO 


. 0012 
. 0025 
. 0040 
. 0080 
. 0115 
0h15 
. 1202 
.0771 


0014 
, 0020 
„ОО 
‚0085 


01.57 


.0581 


‚1155 
.0682 


. 0015 
. 0033 


0047 
0084 
0126 


. 0344 
. 1060 
.0557 


0.0006 
0.0011 
0.0020 
0.0052 
0.0069 
0.0129 
0.0582 
0.1124 


poe 


0.0007 
0.0015 
0.0022 
0.0052 
0.0064 
0.0113 
0.0340 
0.1065 


ee d end 


0.0007 
0.0013 
0.0023 
0.0031 
0.0057 
0.0091 
0.0284. 
0.0985 
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0.7812 
0.3225 
0.1615 
0.1102 
0.0703 
0.01447 
0.0210 


Гооооооо! 


0.8518 
0.3748 
0.1956 
0.1514 
0.0834 
0.0523 
0.0246 


loooooool 


0.8255 
0.3744 
0.197% 
0.1328 
0.0841 
0.0525 
0.0215 


0.8182 
0.57357 
0.1983 
0.1557 
0.0845 
052) 
0215 


[55955551 


8167 
5755 
1987 
1551 
0857 
0522 
. 0211 


[ooooooo] 


IABLE V.- Continued 


AERODYNAMIC-INFLUENCE-CQEFFICIENT MATRICES 


(b) Antisymmetric loadings PN 


Plan form 311 


0.2619 0.1930 
0.7818 0.5867 
1.947} 1.1733 
1.5577 2.4478 
2.1287 5.1126 
1.0690 2.5125 
0.5089 0.9126 
Plan form 512 
0.4566 0.5502 
0.9806 0.7409 
1.9975 1.2558 
1.1026 2.5617 
2.0280 1.7592 
1.0227 2.1767 
0.4856 0.8498 
Plan form 413 
0.4754 O. 3494 
1.0087 0.7750 
2.0136 1.2827 
4.0796 2.3641 
2.0191 1.690 
‚ 1.0168 2.1546 
0.4814 0.8500 


0.4849 
1.0239 


2.0240 


1.0687 
2.0148 
1.0110 
0.4760 


0.4880 
1.0289 
2.0281 
1.0667 
2.0156 
1.0075 
0.4716 


Plan form 511 


0.5601 
0.7912 
1.5009 
2.5691 
1.6690 
2.1587 
0.8201 


Plan form 315 


0. 5680 
0.7974 
1.3076 
2.3720 
11,6616 
2.1511 


0.8231 ` 


0.2248 0.1097 
0.4736 0.2406 
0.8336 0.3969 
1.3018 0.6484 
2.4026 0.9868 
1.8872 1.8992 
1.7910 1.1950 
0.2394 0.1175 
0.4976 0.2534 
0.8564 0.4096 
1.5117 0.6545 
2,5855 0.9817 
1,82 1.8756 
1.7702 1.1576 
0.2485 0 1225 | 
0.5152 0.2620 
0.8727 0.8191 
1.3210 0.6602 
2.576 0.9792 
4. 7865 1.8596 
L. T7517 4.1009 
0.2516 0.1215 
0.5186 0.2651 
0.8787 0.4225 
1.35218 0.6624 
2.357354 0.9181 
1.7716 1.85281 
1.Тіні 1.0865 
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TABLE V.- Continued 


NACA TN 3014 


AERODYNAMIC-INFLUENCE-COEFFICIENT MATRICES 


0.0750 
0.2946 
0.9162 
2.5620 
1.0925 
0.5227 
0.2488 


Plan form 


0.2025 
0.1155 
0.9719 
2.1421 
1.0016 


(b) Continued 


Plan form 321 


0.0612 
0.2156 
0.4963 
1.2307 
2.8688 
1.2702 
0.4821 


222 


0.1568 
0.5185 
0.5598 
1.1552 
2.4972 
1.1052 
0.4171 


0.0907 0.0530 
0.1930 0.0969 
0.3579 0.1645 
0.5925 0.2888 
1.2055 0.721 
2.6290 0.9959 
0.9472 2.5152 


0.1008 0.0479 
0.2108 0.1053 
0.3760 0.1155 
0.5994 0.2917 
1.1795 0.4622 
2.5423 0.9604 
0.9207 2.2511. 
0.1084 0.0515 
0.2246 0.1117 
0.5915 0.1804 
0.6077 0.2948 
1.1629 0.4547 
2.4786 0.9508 
0.8768 2,1721 
0.1118 0.0529 
0.2502 0.111 
0.3977 0.1821 
0.6115 0.2957 
1.1562 0.4504. 
2.4510 0.9155 
0.8566 2.11425 


NACA TN 3014 
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TABLE V.- Continued 


AHRODYNAMTC-INFLUENCE-COEFFICIENT MATRICES 


0.0376 
0.5288 
0.1688 
0.0628 
0.0414 
0.0216 
0.0112 


0.2092 


0.6255 


0.2176 
0.1046 
0.0615 
0.0321 


0.0155 . 


0.2180 
0.6231 
0.2621 
0.1152 
0.0660 
0.0557 
0.0158 


0.2201. 
0.6164 
0.2715 
0.1235 
0.0694 
0.0542 
0.0155 


0.2207 
0.6159 
0.2752 
0.1275 
0.0709 
0.0340 
0.0147 


0.1065 
0.2282 
0.5017 
1.1380 
0.11659 
0.1860 
0.0802 


OOOHROOO 
F 
a 


(b) Continued 


Plan form 334 


0.0657 
0.1509 
0.2568 
0.5480 
1.2783 
0.4856 
0.1549 


Plan form 355 


0.0692 
0.1577 
0.2658 
0.5528 
1.2651 
0.4699 
ο.1ἱμμι 


0.0057 0.0058 
0.0374 0.0191 
0.1129 0.0489 
0.2556 0.1235 
0.7256 0.2545 
1.9242 0.6832 
0.6611 1.8023 


0.0316 0.0142 
0.0661 0.0551 
0.1527 0.0571 
0.2382 0.1129 
0.5908 0.2060 
1.5158 0.5296 
0.5095 1.4097 
0.0364 0.0165 
0.0752 0.0568 
0.1425 0.0608 


0.0805 0.0178 

0.0827 0.0394 

0.1508 0.0635 

0.2551 0.1115 

0.5598 0.1827 

1.221} 0.181 
1 


0.0418 0.0181 
0.0858 0.0100 
0.1542 0.0639 
0.2441 0.1097 
0.5291. 0.175 
1.2915 0.4046 
0.5958 1.1758 
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TABLE V.- Concluded 


AERODYNAMIC-INFLUENCE-COREFFICTENT MATRICES 


OOOOOOO OOOOOOO 


OQOOOOOO 


(b) Concluded 


Plan form 341 


. 0159 
051 
.5381 
1356 
. 0595 
. 025 
.0092 


ооооооо: 


Plan form 342 


Ooooooo 
š 
NO 


-0.0050 


0.0266 
0.1427 
0.8005 
0.2557 
0.0744 
0.0566 


0.0277 
0.0591 
0.1934 
0.6818 
0.2059 
0.0642 
0.0296 


0.0072 -0.0010 


0.0133 
«0526 
‚2488 
.0519 
. 5549 


ооноо 


ОТ 
. 0368 
. 0688 
.2511 
.8162 
.2502 
.0679 


0000000 


Plan form 343 


0.0561 
0.0759 
0.2152 
0.6586 
0.2018 
0.0616 
0.0270 


Plan form 344 


0.0699 
0.2080 
0.5450 
0.1842 
0.0596 
0.0285 
0.0096 


0.0455 
0.0917 
0.2552 
0.6459 
0.2011 
0.0592 
0.0240 


ooo 000 0 
[Ὁ 
η 
n 


0222 


0.0105 
0.0517 
0.0854 
0.5465 
1.2504 


0.5716 


0.0092 
0.0182 
0.0405 
0.0776 
0.2665 
0.9250 
0.2598. 


NACA TN 3014 


0.12758 
0.00042 
0.27772 
0.00050 
0, 37595 
0.00082 
0.87921 
-0.03954 


TABLE VI 


1 -1 
VALUES OF [eos nd] | {sin négp FOR INFLUENCE-COEFFICIENT CALCULATIONS 


(a) Symmetrical distributions 


[гов "UB { sin no, b for - 


0.05490 
0.15587 
0.12089 
0.19161 
0.17232 
0.35209 
0.67744. 


-0.52044 


0.10825 
0.07457 
0.24856 
0.10198 


0.45929 
0.51266 


-0.60008 
~0.04983 


0.00112 
0.36868 
0.00289 
0.72229 
0.0806), 
-0.54349 
-0.00127 


-0. 01658 


πτος NL VOVN 


£9 


TABLE VI.- Concluded 
4-1 
VALUES OF | [eos n5] | (ein пӨо» FOR INFLUENCE-COEFFICIENT CALCULATIONS 


(b) Antisymmetrical distributions 


-l 
| [eos ni] | { sin no for - 


ast eoa оаза teet = oet eoo ret 0-б]; eo et otn = 0.9 


0.02547 
0.00008 
0.06011 
0.00013 
0.10633 
0.00029 


0.45950 


0.01617 
0.04767 
0.05925 
0.06915 
0.07311 
0.19013 
0.53108 


0.05413 
0.03801 


0.13453 
0.06132 


0.13699 
0.06636 


-0.04053 
-0.61138 


ητος NL VOVN 
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.19509 
. 58268 
22221 
ТОТ 
8314-7 
. 92388 
.98079 
. 00000 


. 01877 
. 13889 


20787 
24520 


.21520 


20787 
15889 
0#877 


00000 
98079 
92588 
83147 
70711 


«22251 
. 58268 


„19209 


.06250 


TABLE VII 


MATRICES USED IN LIFI-DISTRIBUTION CALCULATIONS 


0.22221 


0.12500 


(a) Symmetrical distributions 


[sin nə, | = 


0.83147 
0.92388 
0.19509 
~O. 70711 
-0.98079 
-0. 38268 
0.55557 
1.00000 


0.13889 
0.24520 
0.04877 
-0.20787 
-0.20787 
0.01877 
0.24520 
0.13889 


5.00000 
2.77785 
-1.91340 
-4.90395 
-5.53555 
0.97545 
4.61940 
4.15735 


0.12500 


0.98079 0.98079 
0.38268 -0.38268 
-0.83147 -0.83147 
-0.TOT11 0.70711 
0.525551 0.55557 
0.92588 -0.92588 
-0.19509 -0.19509 
-1.00000 1.00000 
[515 neg] = 
0.17678 0.20787 
0.17678 0.04877 
-0.17678 -0.24520 
-0.17678 0.13889 
0.17678 0.13889 
0.17678 -0.24520 
-0.17678 0.04877 
-0.17678 0.20787 


| cos non | In] = 


7.00000 
1.56565 
-6.46716 
-3 . 88899 
1.91977 
5.82029 
-2.67876 
-6.86555 


9.00000 
-L. 75581. 
-8.31492 

5.00015 

6.56599 
-7.48323 
-3.412 

8.82711 


ΕἾΝ 


0.12500 


0.12500 


0.85187 
-0.92388 
0.19509 
0.70711 
-0.98079 
0.38268 
0.55557 
-1.00000 


0.23097 
-0.09567 
-0.09567 

0.25097 
-0.23097 

0.09567 

0.09567 
-0.23097 


11.00000 
-6.11127 
-11. 20948 
10.78869 
-7.77821 
-2.14599 
10.16268 
-9.14617 


0.12500 


0.55557 
-0.92588 
0.98079 
-0.70711 
0.19509 
0.58268 
-0.83147 
1.00000 


0.24520 
-0.20787 
0.13889 
-0.04877 
-0.04877 
0.13889 
-0.20787 
0.24520 


15.00000 


-10.80911 


1.978 
2.55617 
-9.192145 
12.75027 


-12.010)-4 


T.22211 


0.12500 
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TABLE VII.- Continued 
MATRICES USED IN LIFT-DISTRIBUTION CALCULATIONS 
(a) Concluded 
gin no, ЫШ 
EE Ше 
1.00000 8.51515 21.30990 35.19159 5.21655 16.88178 37.02025 15.00000 
| 1.00000 7.24272 12.07120 7.00000 -9.00000 -26.55640 -31.38512 -15.00000 
1.00000 5.29620 1.75580 -10.47641 -13.46967 5.86276 22.95020 15.00000 
1.00000 2.00000 -5.00000 -7.00000 9.00000 11.00000 -13.00000 -15.00000 
1.00000 0.70589 -5.89795 1.67719. 6.01555 -12.97549 2.05019 15.00000 
1.00000 -1.24263 -2.07105 7.00000 -9.00000 4t ,55651. 5.38473 -15.00000 
1.00000 -2.54328 2.85220 -1.39237 -1.79019 6.23084 -11.02088 15.00000 
1.00000 -3.00000 5.00000 -7.00000 9.00000 -11.00000 13.00000 -15.00000 
[às] - 
5,1258 -1.8481 О -0.1514 -0.0481 0 -0.0165 
-0.9h22 2.6151 -1.0195 0 -0.1017 ο -0.0411 0 
0 -0.7022 1.7999 -0.7191 O -0.0815 0 | -0.0226 
-0.0417 0 -0.56h9 1.4140 -0.5739 o -0.0715 ο 
ο -0.0468 0 -0.4881 1.20217 -0. 499% . О -0.0506 
-0.0102 О -0.0490 O -0.4494 1.0824 -0.4814 0 
O -0.0160 O -0.0557 0 -0.4555 1.0196 -0.1106 
-0.0065 O -0.0251 ο -0.0842 O -0.8055 1.0000 
[5] - 

0.059276 -0.018861 0.011695 -0.007815 0.005220 -0.003236 0.001554 
-0.018861 0.050969 -0.02667 0.016915 -0.011048 0.006775 -0.005256 
-0.027584 -0.007813 0.044496 -0.022097 0.013246 -0.007813 0.003666 

0.011048 -0.034057 -0.003236 0.040830 -0.018861 0.010158 -0.004577 
-0.006473 0.015625 -0.057723 0 0.037723 -0.015625 0.006475 

0.004577 -0.010138 0.018860 -0.040830 0.003236 0.034056 -0.011048 
-0.003666 0.007815 -0.013246 0.022097 -0.044496 0.007813 0.027584 

0.005256 -0.006775 0.011048 -0.016913 0.026674 -O 050969 0.018861 
-0.00155h 0.002236 -0.005220 0.007813 -0.011693 0 018861 -0.059276 
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TABLE VII.- Continued 
MATRICES USED IN LIFT-DISTRIBUTION CALCULATIONS 
(b) Antisymmetrical distributions 
[sin n6, | = 
0.58268 0.70711 0.92388 1.00000 0.92588 
0.70711 1.00000 0.70711 0.00000 -0.70111 
0.92588 0.70711 -0.38268 -1.00000 -0.38268 
1.00000 0.00000 -1.00000 0.00000 1.00000 
0.92588 -0.70711 -0.38268 1.00000 -0.38268 
0.70711. .-1.00000 0.70711 0.00000 -0.701711 
0.58268 -ο.τοτιι 0.92588 -1.00000 0.92588 
[ein non] шы 
0.09567 0.17678 0.25097 0.25000 0.25097 
0.17678 0.25000 0.17678 0.00000 -0.17678 
0.25097 0.17678 -0.09567 -0.25000 -0.09567 
0.25000 0.00000 -0.25000 0.00000 0.25000 
0.25097 -0.17678 -0.09567 0.25000 -0.09567 
0.17678 -0.25000 0.17678 0.00000 -0.17678 
0.09567 -0.17678 0.25097 -0.25000 0.25097 
| сов no, | In| = 
2.00000 4.00000 6.00000 8.00000 10.00000 
1.81776 2.82844 2.29608 0.00000 -3.82680 
1.3122 0.00000 -h.2h266 -8.00000 -~-7.07110 
0.76556 -2.82844 -5.54348 0.00000 9.23880 
0.00000 -і.00000 0.00000 8.00000 0.00000 
-0.76536 -2.82844 5.54328 0.00000 -9.25880 
-1.41422 0.00000 4.24266 -8.00000 7.07110 
-1.84776 2.828: -2.29608 0.00000 5.82680 
-2.00000 4.00000 -6.00000 8.00000 -10.00000 
ln] - 
0.06250 | 
12500 
0.12500 
0.1250 
0.12500 
0.12500 


0.70711 0.38268 
-1.00000 -0.70711 
0.7071 0.92588 
0.00000 -1.00000 
-0.70711 0.92588 
1.00000 -0.70712 
-ο. 70711 0.38268 
0.17678 0.09567 
-0.25000 -0.17678 
0.17678 0.23097 
0.00000 -0.25000 
-0.17678 0.23097 
0.25000 -0.17678 
-0.17678 0.09567 
12.00000 1}.00000 
-8.48532 -12.93432 
0.00000 9.89954 
8.18552 -5.59722 
-12.00000 0.00000 
8.18552 2.22122 
0.00000 -9.89954 
-8.48532 12.93432 
12.00000 -1}.00000 
0.12500 
0.12500 
0.06250 
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.050056 
.020415..0. 
.025981. 
.009566 -0. 
.005666 ο. 
.002595 
001450 0. 
. 000644 


5.92512 
5.69556 
5.52594 
2.82842 
2.22228 
1.53074 
0.78036 


1 
OOOOOOW! 


.020#15 
054108 
-011048 
050649 
011959 
.006097 
003236 
. 001430 


TABLE VII.- Concluded 


MATRICES USED IN LIFT-DISTRIBUTION CALCULATIONS 


1.49812 
10.45260 
5.09116 
-0.00000 
-3.40172 
-4.32956 
-2 . 88384 


0.0000 
-0.0065 


0.018062 
-0.031464 
0.049442 
-0.008456 
-0.032081 
0.012605 
-0.006097 
0.002595 


28.113396 
11.08668 
-.13292 
-8.18526 
-2.76150 

%,59222 

5.65188 


(о) Concluded 


E = uj - 


sin 9 


41.00672 
0.00009 
-14.39992 
0.00000 
9.62152 
0.00000 
-8.15672 


[s] - 


-0.0110A8 
0.025459 
-0.028872 
0.048012 
-0.007815 
-0.022081 
0.011959 
-0.004666 


. 009396 
.019505 
.022028 
.028228 
. 018012 
.008156 
030649 
. 009366 


7.55660 
-18.47780 
-6.88820 
14.1h210 
_-8..60250 
-7.65370 
9.43980 


0.0000 
-0.0898 
0.0000 
-0.5576 
1.2027 


20,4199 


0.0000 
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hlz.h9h36 27.46184 
-31.35780 -25.86892 
15.27348 23.28158 

0.00000 -19.79894 
-10.20516 15.55596 

12.98868 -10.71518 

-8.65152 5.46252 

-0.0226 0.0000 

0.0000 -0.0166 
-0.0619 0.0000 
0.0000 -0.0596 
-0.4666 0.0000 
1,0824 -0.5896 

-0.3670 1.0196 
-008456 0.007966 -O 
.017562 -0.016269 0 
.018861. 0.017362 -0 
022028 -0.019505 ο 
.028871 0.023459 -0 
.049442 -0.031464 ο 
.011048 0.054108 -O 
025984 -0.020415 ο 


y* 
т 
y* - af 
i/(y* - ww) 
l/(y* + q*) 
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TABLE VLII.- Continued 
' COMPUTING PROCEDURE FOR F FUNCTIONS 


(b) Bow 2 


0.00579 0.0032} 0 0.15565 
5.70.51 5.62777 5.81422 2.18877 
0.85555 0.85555 0,85555 0.85555 


at = ух, (22) = tan А9). 
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0.19728 
2.41980 
0.30866 


7 |®® ®® @® ® ® ® ®® OOO ®® © OO @©) © Ὁ 


enm) rm. 


(e) Row 5 
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7.59 
5.4247 5. 
5,5866 
2.0610 
1.2464 
0.7125 
0.4229 


TABLE IX 


SYMMETRIC AND ANTISYMMETRIC P MATRICES FOR PLAN FORM 255 


-9.0860 
-10.2951 
-10.9636 
-10.7385 
-9.9951 
-9.0643 
-8.1161 
-7.2196 


8.3357 
9650 
2.6601 
2.1859 
1.5057 
0.7409 
0.5515 


-8.3357 
-9.7155 


-10.6795 
-10.6015 


-9.9214 
-9.0182 
-8.0817 
-7.1891 


10.4164 
1.8390 
1.6921 
2.6370 
1.5099 
0.8560: 
0.5726 


12.6828 
10.6672 
7.1238 
5.7162 
1.9676 
1.0551 
0.4498 


al F | = 


14.1611 14.9907 
12.709} 15.7911. 
10.175 11.9020 
6.2896 9.2061 
5.0125 5.5952 
1.31254 2.3579 
0.5975 0.8970 
[τα] = 
-5.0097 -2.3292 
-1.9575 -2.854% 
-7.1258 -5.9851. 
-8.9562 -6.2896 
-9.1347 -7.9649 
-8.5679 -8.0014 
-7.7636 -7.4021 
-6.9106 -6.6056 


15.44.55 
14.5488 
12.7114 
10.6338 
7.9725 
44.096 
1.5789 


15.6705 
14.6188 
13.0781 
11.2167 
9.0854 
6. T2} 
3.0437 


-0.8120 
-0.9156 
-1.1327 
-1.5732 
-2.5519 
~ 4096 
-5.5165 
-5.1811 
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TABLE X 


INTEGRATING MATRICES FOR LOAD AND MOMENT COEFFICIENTS 


kok 


|0.01915 0.05757 0.05454 0.06942 0.08165 0.09070 0.09629 0.04909 | 


7 


[0.01878 0.03472 0.04555 0.04915 0.04528 0.03487 0.01833 0.00196 | 


Pal 


[0.01995 0.05596 0.05715 0.06556 0.08156 0.08164 0.11565| 


[ο 


|0.009% 0.01755 0.02267 0.02454 0.02267 0.01736 0.00939 | 


0.01035 
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(a) Symmetrical lift distributions. 
(b) Antisymmetrical lift distributions. 


Figure 1.- Spanwise lift distributions for wings of very low aspect 
ratio (from ref. 3). 
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Spanwise loading coefficlent, 


(c) Lift distribution for inboard flap. 
(а) Lift distribution for outboard aileron. 


Figure 1.- Concluded. 
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(a) Symmetrical lift distributions. 


(b) Antisymmetrical lift distributions. 


Figure 2.- Spanwise lift n for plan form 511 (A = 1.5; 
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(c) Lift distribution for inboard flap. 
(4) Lift distribution for outboard aileron. 


Figure 2.- Concluded. 
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(a) Symmetrical lift distributions. 


(b) Antisymmetrical lift distributions. 


Figure 3.- Spanwise lift distributions for plan form 312 (A = 1.5; 
À = 0.25). 
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(c) Lift distribution for inboard flap. 
(а) Lift distribution for outboard aileron. 


Figure 4.- Concluded. 
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(a) Symmetrical lift distributions. 
(b) Antisymmetrical lift distributions. 


Figure h.- Spanwise lift distributions for plan form 313 (A = 1.5; 
À = 0.50). 
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(с) Iift distribution for inboard flap. 
(d) Lift distribution for outboard aileron. 


Figure 4.- Concluded. 
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(a) Symmetrical lift distributions. 


(b) Antisymmetrical lift distributions. 
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Figure 5.- Spanwise lift distributions for plan form 314 (A = 1.5; 
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(c) Lift distribution for inboard flap. 


(d) Tift distribution for outboard aileron. 


Figure 5.- Concluded. 
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(a) Symmetrical lift distributions. 
(b) Antisymmetrical lift distributions. 


Figure 6.- Spanwise lift distributions for plan form 315 (A = 1.5; 
^ = 1.50). 
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(c) Lift distribution for inboard flap. 


(d) Lift distribution for outboard aileron. 


Figure 6.- Concluded. 
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(a) Symmetrical lift distributions. 
(b) Antisymmetrical lift distributions. 


Figure 7.- Spanwise lift distributions for plan form 321 (A = 3.0; 
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(c) Lift distribution for inboard flap. 


(d) Lift distribution for outboard aileron. 


Figure 7.- Concluded. 
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(a) Symmetrical lift distributions. 
(b) Antisymmetrical lift distributions. 


Figure 8.- Spanwise lift distributions for plan form 522 (A = 5.0; 
À = 0.25). 
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(ο) Lift distribution for inboard flap. 
(d) Lift distribution for outboard aileron. 


Figure 8.- Concluded. 
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Figure 9.- Spanwise lift distributions for plan form 323 (A = 3.0; 
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(а) Symmetrical lift distributions. 


(b) Antisymmetrical lift distributions. 


= 0.50) 


97 


NACA TN 3014 


i3A 


cote 
ELL WAA ДИ 
AENEA 
πα AMIE 


-HMUCHE 
Е ИДУ ΠΗ} 


SN ETATE 
Е ГИДА АЕ 
Рур: 
Е ИНЕТ У 
LL y: 
TRIN NUMINE 


«e 
i -і 


Ὅτ 
: 29 '*uorogreoo Зитреот espwuredg 
29 


PASSES 


анала ЫН 
TITAN CME 


dud 
SEAL 


- 


© 
- 


*querogieoo Зитреот эвтицзая 


P za spanwise ordinate, 2 


(ο) Lift distribution for inboard flap. 


(4) Lift distribution for outboard aileron. 


Figure. 9.- Concluded. 
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(&) Symmetrical lift distributions. 
(b) Antisymmetrical lift distributions. 


Figure 10.- Spanwise lift distributions for plan form 52h (A = 5.0; 
= 1.00). 
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(c) Lift distribution for inboard flap. 
(a) bas distribution for outboard aileron. 


Figure 10.- Concluded. 
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Figure 11.- Spanwise lift distributions for plan form 225 (A = 5.0; | 
о). 
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(а) Symmetrical lift distributions. 


(о) Antisymmetrical lift distributions. 
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(с) Lift distribution for inbosrd flap. 


(d) Lift distribution for outboard aileron. 


Figure 11.- Concluded. 
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(a) Symmetrical lift distributions. 


(b) Antisymmetrical lift distributions. 


Figure 12.- Spanwise lift ΕΓ plan form 331 (A - 6.0; | 
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Dimensionless spanwise ordinate, y* 
(c) Lift distribution for inboard flap. 


(а) Lift distribution for outboard aileron. 


Figure 12.- Concluded. 
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(&) Symmetrical lift distributions. 
(b) Antisymmetrical lift distributions. 


Figure 15.- Spanwise lift distributions for plan form 332 (A = 6.0; 
Aog Оо 
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(ο) Idft distribution for inboard flap. 
(d) Lift distribution for outboard aileron. 


Figure 15.- Concluded. 
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(&) Symmetrical lift distributions. 
(b) Antisymmetrical lift distributions. 


Figure 14.- Spanwise lift distributions for plan form 333 (A = 6.0; 
À = 0.50). EM | 
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Dimensionless spanwise ordinate, y* 


(c) Lift distribution for inboard flap. 
(а) Lift distribution for outboard aileron. 


Figure 1+.- Concluded. 
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(a) Symmetrical lift distributions. 
(b) Antisymmetrical lift distributions. 


Figure 15.- Spanwise lift distributions for plan form 334 (A = 6.0; 
A= 1.00 
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(c) lift distribution for inboard flap. 
(а) Lift distribution for outboard aileron. 


Figure 15.- Concluded. 


109 


110 NACA TN 3014 


в ETTI 5 
: ЕЕЕ ЕЕ 
ЖИК ИШ ШИИ ШШ T Б ШО Ee 

8 SEEEREEEEEEREEEEE EE 
: ШЕ [η 

4 


„= ΗΝ 
Е ΠΗΓΗ 
ΠΠ 


р Е 

ЕРМЕН Е В ТЕ 
ЕЕ KILI: 
LIT TILL TIN 
NNNM Е T ж аны ШЕ 


[UP 


Spanwise loading coefficient, 


а εἶνε. eae os a 


(a) Symmetrical lift distributions. 
(b) Antisymmetrical lift distributions. 


Figure 16.- Spanwise lift distributions for plan form 335 (A = 6.0; 
À = 1.50) 
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(c) Lift distribution for inboard flap. 


(а) Lift distribution for outboard aileron. 


Figure 16.- Concluded. 
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(&) Symmetrical lift distributions. 
(b) Antisymmetrical lift distributions. 


Figure 17.- Spanwise lift m for plan form 341 (A = 12.0; 
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(c) Lift distribution for inboard flap. 
(à) Lift distribution for outboard aileron. 


Figure 17.- Concluded. 
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(a) Symmetrical lift distributions. 
(b) Antisymmetrical lift distributions. 


Figure 18.- Spanwise lift distributions for plan form 342 (A = 12.0; 
À = 0.25). 
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(c) Lift distribution for inboard flap. 
(а) Lift distribution for outboard aileron. 


Figure 18.- Concluded. 
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(a) Symmetrical lift distributions. 
(b) Antisymmetrical lift. distributions. 


Figure 19.- Spanwise lift distributions for plan form 555 (А = 12.0; 
^ = 0.50) 
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Spanwise loading coefficient, 


(c) Lift distribution for inboard flap. 
(а) Lift distribution for outboard aileron. 


Figure 19.- Concluded. 
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(a) Symmetrical lift distributions. 
(b) Antisymmetrical lift distributions. 


Figure 20.- Spanwise lift distributions for plan form 344 (A = 12.0; 
- 1.00) 
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Dimensionless spanwise ordinate, y* 
(c) Lift distribution for inboard flap. 
(d) Iift distribution for outboard aileron. 


Figure 20.- Concluded. 
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Effective angle of attack for unit flap deflection, сср 
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Figure 21.- Variation of flap effectiveness with flap-chord ratio. 
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